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ABSTRACT 


A computer model of an aviation gas turbine test 
facility was developed and analyzed to compare computer 
generated gas pressures, temperatures, and velocities to 
those of an actual, operational test facility. The 
computer model simulated the Royal Danish Air Force Test 
Facility at the Skrydstrup Air Force Base, Denmark. The 
Parabolic, Hyperbolic, or Elliptic Numerical Integration 
Code Series (PHOENICS) was used. The Computer Model output 
compared well with actual data obtained from the Denmark 


test facility. Tabular Data and Figures are presented. 
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I. INTRODUCTION 


A. BACKGROUND 

The United States Navy desires to provide standardized 
test facilities for the testing of its aviation gas turbine 
jet engines (Ref 1). 

Today’s high performance and high technology jet 
engines require accurate, consistent, and controlled 
testing to assess performance throughout the jet engine’s 
entire operational range Modern day test facilities must 
incorporate state of the art equipment and instrumentation 
such that jet engines may be tested in or out of the 


aircraft, around the clock, and year round. 


B. NATURE OF THE PROBLEM 

Presently, 75 Navy and Marine Corp activities worldwide 
operate 81 Jet Engine Test Facilities (Ref 1). Of these, 
60 are fully operational and 51 are over 25 years old. 
Immediate problem areas which hinder operational needs 
include: structural degradation, obsolescence, noise, and 
air pollution. 

1. Structural Degradation 

High jet engine exhaust gas temperatures and water 

impingement on diffusers and interior concrete walls have 


caused cracking and spalling of test facility interiors. 


Concrete reinforcing rods exposed to a hot, moist 
environment will corrode and thereby reduce structural 
integrity. Noise reduction baffles and exhaust gas turning 
vanes are being damaged by flying concrete debris. 
Replacement of baffling and turning vanes is required as 
often as every three years. The maintenance requirements 
navy wide has significantly increased the workload and 
operational costs of the Navy Public Works Department. 
2. Obsolescence 

Several of the test facilities require major 
modernization or even complete replacement due to current 
aircraft designs and test requirements. Examples of 
modernization include improved engine handling systems and 
engine thrust stands. Increased access to engine 
components and test points greatly improves testing 
efficiency. Additionally, with today’s high technology 
engines, state of the art data acquisition systems are 
required. 

3. Noise 

Most of the older test facilities use water for 
exhaust gas cooling. The water and the high exhaust gas 
temperature deteriorates the acoutic baffling silencing 
capabilities. The increased . noise level degrades the work 
environment. Worker efficiency is reduced and the high 
test facility interior sound levels are being suspected of 


contributing to engine deterioration. State of the art 
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engines using thinner, lighter materials are more sensitive 
to acoustic fatigue (Ref 1). 

a. hale Pollution 

Today’s high power and high mass flow engines 

constitute a major problem to the environment. Visible 
particulate matter is emitted in excess of allowable 
environmental pollution standards. By Executive Order 
Number 11282, May 26, 1966 (Ref 2), Jet Engine Test 
Facilities must meet local laws concerning the environment. 
Judicial action by the state of California against Jet 
Engine Test Facilities within California (Ref 3) has 
identified a need to address pollution control measures 
during the design, construction, and operation of Jet 
Engine Test Facilities. 
eee, LYPES OF TEST FACILITIES 

The Navy operates five types of Jet Engine Test 
Facilities (Ref. 1). Each facility utilizes some form of 
hold down device for the engine or the aircraft and 
instrumentation to monitor the performance of the engine. 
Each facility must be designed such that several engine 
types may be tested (e.g., turboprop, turbofan, turboshaft, 
mimeboj7et). The types of facilities are: 

fi) Run-up Pad 

The simplest test facility is the run-up (power 

check) pad. Normally, the engine is tested while installed 


in the aircraft. Therefore, the aircraft acts as the hold- 


CA 


down device. The run-up pad is usually located in remote 
locations due to high noise levels encountered during 
testing. 
2. sound Suppressors 
These are usually equipment packages which are used 
to reduce inlet and exhaust noise from a tied down 
aircraft. The sound suppressors may be used to augment or 
to substitute for the run-up pads. 
3. Open Test Stand 
This is an outdoor test facility which supports and 
constrains the engine during test. Due to noise levels, 
the open test stand must be located in a remote area. 
Portable sensors and instrumentation are utilized to 
monitor engine performance. 
4. Hush House 
This is a total aircraft enclosure. It permits 
round the clock and year round testing. Normally on twin 
engine aircraft, only one engine may be tested in full 
after burner. The other engine is held at idle. 
5. Fully Enclosed Test Cell 
This test facility is normally capable of testing 
several types and sizes of jet engines. The engines may be 
in or out of the aircraft. The facility is normally 
equipped with all the necessary instrumentation and data 


gathering equipment required to perform routine or complex 


testing after jet engine rework or overhaul. Testing may 
be performed round the clock and year round. 

The type of test facility to utilize at an activity is 
dependent upon several factors. Availability of military 
construction funds, required control of noise and emissions 
(environmental concern), amount of maintenance testing, and 
military and civilian manpower requirements are all key 


factors. 


D. RESEARCH GOAL 

To provide quality jet engine testing, to comply with 
local environmental standards, and to ensure fleet 
readiness, the Navy embarked on a comprehensive program in 
1984 to upgrade and to modernize as necessary all jet 
engine test facilities under its cognizance. To this end, 
it was necessary to review Seen jet engine test 
procedures and to examine jet engine test facilities to 
determine which facilities met desired criteria. 
Additionally, the Navy decided to explore the feasibility 
of the computer modeling of aviation gas turbine jet engine 
test facilities. It was felt that if a computer model 
could accurately represent an actual operational test 
facility then a computer model could also be used to 
predict the operation of future proposed test facilities. 
The purpose of this research was to computer model the 


Royal Danish Test Facility and to compare computer 


generated temperatures, pressures, and velocities with 
actual data obtained in the field and to investigate 
techniques necessary to make computer simulations of flow 


and temperature fields useful for predictive purposes. 
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II. TEST FACILITY OPERATION 


The basic geometry of the Denmark test facility is 
shown in Figure C.1. For simplicity all storage areas, 
engine preparation rooms, and control modules are not 
included. The test facility is basically comprised of the 
primary air inlet, a test bay, an exhaust duct (also called 
an augmentor tube), secondary air inlets, and an exhaust 
stack. Turning vanes, which are not shown, are located in 
the exhaust duct to exhaust stack transition region. This 
feature better utilizes the total exhaust stack cross 
section with an attendant reduction in exhaust gas noise 
(Ref 4). Baffles (not shown) are located at three places: 
At the top of the exhaust stack, within the secondary air 
inlets, and within the primary air inlet (test bay doors). 
The baffles also provide noise reduction (Ref 4). 

During a jet engine test, the engine is aligned with 
the centerline of the exhaust duct. This can be 
accomplished either with the engine in the aircraft or with 
the engine supported in a test stand. With proper 
alignment, a near uniform flow profile develops (Ref 3) 
inside the augmentor tube. Primary air iS drawn in through 
the primary air inlet (test bay door/baffle configuration). 
For a F-100 Engine, at full throttle and after burner on, 


the primary air mass flow rate is 260 kg/sec of which 
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100 kg/sec 1s the mass flow rate through the engine 

(Ref. 4). As the exhaust gas enters the augmentor tube, 
air from within the test bay is entrained into the 
augmentor tube. The exhaust gas moving through the 
augmentor tube also entrains air coming through the 
secondary air inlets. The secondary air mass flow rate is 
340 kg/sec (Ref. 4). Both the entrained primary and 
secondary air play an important role in the test facility 
operation. Hot and noisy jet engine gases must be treated 
prior to leaving the test facility in order to comply with 
environmental standards. To treat jet engine exhaust, its 
temperature, noise level, and particulate density must be 
lowered. 

The ratio of total "cooling" air (primary and 
secondary) mass flow rate to the engine mass flow rate is 
called the Augmentation Ratio. The Augmentation Ratio is a 
key to successful test facility operation. With 
appropriate primary air and secondary air mass flow rates, 
the exhaust gas is sufficiently cooled which helps prolong 
duct/stack life and the exhaust gas is sufficiently diluted 
such that proper pollution levels are obtained (Ref. 2). 
The Augmentation Ratio of the Denmark Test Facility is 5:1. 

With a computer model that can accurately predict test 
facility operation, design personnel will have the ability 
to easily change model parameters such as engine location 


or structural dimensions and to observe their effect on 


test facility operation. For example by changing the 
secondary air inlet configuration, the augmentation ratio 
is changed. This in turn has an effect over exhaust gas 
cooling. Additionally, design engineers may wish to 
experiment with various exhaust stack heights or exhaust 
duct openings. A mere changing of a few lines of computer 
model source code accomplishes this. In all, a computer 
model could save significant time and money in the design 


of future jet engine test facilities. 


Iitl. DESCRIPTION OF PHOERNMeS 


A. WHAT IS PHOENICS 

PHOENICS is a general-purpose, computer-code system for 
Simulating fluid flow, heat transfer, and chemical reaction 
processes in engineering equipment or the environment. 
PHOENICS is an acronym for Parabolic, Hyperbolic, or 
Elliptic Numerical Integration Code Series. Typical 
applications of PHOENICS include: Thermal hydraulics of 
nuclear reactors and other power plants; external 
hydrodynamics and aerodynamics of cars, ships, buildings, 
etc; air movement within an enclosed space, both under 
normal conditions and in the presence of a fire; flows in 
pumps and turbomachinery. PHOENICS was developed by D. 
Brian Spalding of CHAM Ltd and of Imperial College of 
London, England. Information regarding the history and 
development of PHOENICS is contained in the user’s guide 


(Ref 5). 


B. HOW PHOENICS WORKS 

The PHOENICS computer code consists of two separate 
computer programs: EARTH and Ql. The EARTH program 
resides in the computer mainframe as an object code and 
cannot be accessed by the user. The user couples his 


problem (defined by the user) to the EARTH code. The EARTH 
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program contains the main flow simulating software and 
various coding sequences which simulate the Laws of Physics 
aS applied to the conservation principles: energy, mass, 
momentum. The Ql program defines the problem and differs 
from one flow simulation program to the next. The Qi 
program represents the user’s decision as to what problem 
is to be simulated. The Q1 program activates the EARTH 
program which solves the problem. 

The PHOENICS code is capable of describing flow 
Simulations in terms of distributions in time and space. 
The distributions involve assigning numerical values for 
all dependent variables to an orderly array of cells 
located in the physical domain of the problen. 
Temperature, pressure, concentration, KE (Turbulence 
Kinetic Energy), and EP (Kinetic Energy Dissipation Rate) 
values are located at a cell’s center while velocity values 
are located at the positive faces of the cells. A cell 
configuration in the X-Y plane is illustrated in Figure 
C.2. Temperatures, pressures, and concentrations would be 
located at N, S, E, W, and G for example while velocities 
would be at N’, S’, E’, and W’. A similar procedure is 
utilized for the Z (low-to-high) direction. This 
convention is called the "STAGGERED GRID" arrangement. If 
the flow simulation is time dependent, the distributions 


are obtained for each variable at a succession of times 


il 


instants. The Denmark test facility model flow simulation 
is for steady state conditions. 

In general, a PHOENICS flow simulation is four 
dimensional: 1 time and three space coordinates. For ease 
of reference to the user, PHOENICS measures time in the 
past-to-future direction; X coordinates in the west-to-east 
direction, Y coordinates in the south-to-north direction, 
and Z coordinates in the low-to-high direction. This 
concept is illustrated in Figure C.3. Although Figure C.3 
represents a Cartesian configuration, individual cells may 
vary in shape and size as determined by a user defined 
code. This option, called BFC, was not utilized in this 
investigation. The Denmark test facility model is a 
Cartesian configuration of cells which vary in size. 

PHOENICS calculates the dependant variables over the 
complete set of cells, which defines the operational (e.g., 
system, physical) domain, with linear algebraic equations 
that represent the conservation laws. Here the partial 
differential equations are those governing the conservation 
of energy, momentum, and mass and the linear algebraic 
equations are obtained by the finite volume method. The 
finite volume method is a method having some features from 
the finite differencing method and the finite element 
method. The conservation partial differential equations 
will not be discussed here; except to say that the original 


non linear, coupled conservation partial differential 


(B2: 


equations are linearized. Interested readers are directed 
to the PHOENICS User’s Manual (Ref. 5) for detailed 
Piimearmnation. 

PHOENICS utilizes an iterative process involving a 
multi-stage sequence of dependent variable value 
adjustments. The process is repeated several times until 
the left and right imbalance (called a residual) of each 
equation is reduced to a negligible amount. A variety of 
iterative procedures are available to the user. The 
"SLAB-SWEEP" iterative process feature of PHOENICS for 
equation imbalance reduction was chosen for the Denmark 
test facility model. Slabs are arrays of cells having the 
same Z (low-to-high) coordinate. PHOENICS conducts many 
mathematical operations over a slab in an attempt to reduce 
dependent variable imbalances among cells within the slab 
prior to moving on to the next slab. The adjustment 
operations are referred to as "SLABWISE SOLUTIONS". The 
operations begin at the lowest numbered slab and continue 
sequentially to the highest numbered slab. A sweep is a 
completed set of slabwise operations. The equations for 
dependent variable values at one slab are coupled to values 
of other dependent variables at neighbor slabs. For this 
reason numerous sweeps are normally required until all 
dependent variable imbalances are negligible. This is a 
necessary but not a sufficient condition for problem 


convergence. 


13 


GC. PHOENICS INPUT LANGUAGE 

The source code used in a PHOENICS Q1 program is 
arranged in a 24 group data structure. Figure C.4 lists 
the group categories of the data structure. The developers 
of the PHOENICS code suggest this ordering of the groups 
normally conforms to a well thought out flow simulation 
problem description. However, any group ordering the user 
desires will be acceptable to PHOENICS. It is the PHOENICS 
defined variables and keywords plus the user defined 
variables that actually trigger desired PHOENICS responses 
aS opposed to a particular group ordering. A detailed 
description of group functions is explained in the PHOENICS 


User’s Manual (Ref. 5). 


D. MODEL GEOMETRY | 

PHOENICS utilizes three distinct procedures for 
defining a flow simulation operational domain (grid): 
Cylindrical-polar, Body-fitted, and Cartesian coordinate 
systems. 

1. Cylindrical-Polar Coordinate System 

The cylindrical-polar grid arrangement consists of 

planes with constant Z coordinates that are perpendicular 
to an axis of rotation, planes with a constant X coordinate 
(usually called the Theta coordinate) that intersect with 


the axis of rotation giving the X coordinate a 
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representation in radians, and concentric cylindrical 
surfaces with a constant Y radial coordinate. 
2. Body-Fitted Coordinate System 

This arrangement is best thought of as starting 
with an orthogonal group of cells made of clay. The cells 
are then squeezed, bent, twisted, and stretched into the 
desired shape for the fluid flow simulation operational 
domain. Each cell although distorted in shape and size 
still remains in contact with its original neighbors. 

3. Cartesian Coordinate System 

The Cartesian grid arrangement is composed of cells 
in a mutually perpendicular X, Y and Z coordinate system. 
As such the cell faces form planes that intersect at right 
angles. The spacing between the cell faces (planes) is 
user defined. 

Due to the regularity of a Cartesian coordinate 
system, the information required to specify the geometry is 
rather modest. The PHOENICS variables NX, NY, and NZ 
specify the number of cells the operational domain is to 
have in the X, Y, and Z directions respectively. The 
PHOENICS variables XULAST, YVLAST, and ZWLAST specify the 
domain lengths in the X, Y, and Z directions. The "METHOD- 
OF-PAIRS" is a PHOENICS coding sequence to specify the 
number of cells per length of domain in the X, Y, and Z 
Girections. XFRAC, YFRAC, and ZFRAC are PHOENICS variables 


for arrays used in the method-of-pairs coding sequence. 


= 
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The odd indexed array elements, for example, XFRAC(1), 
XFRAC(3), XFRAC(5), etc., specify the number of cells ina 
particular length of operational domain while the even 
indexed array elements XFRAC(2), XFRAC(4), and XFRAC(6) 
specifies the length of each cell. A negative argument for 
the first element in each array tells PHOENICS that the 
method-of-pairs is being utilized. For example, "XFRAC(1) 
=-2.0;XFRAC(2)=0.1" tells PHOENICS the method-of-pairs is 
to be used and that the first cell group in the X direction 
is to have two cells and the length of each cell is 0.1 
meter. "XFRAC(9)=4.0;XFRAC(10)=0.15" specifies the number 
of cells and lengths of each cell in the "FIFTH" xX 


direction cell group. 


E. POROSITY AND BLOCKAGE 

In any flow simulation problem there exists walls, 
submerged bodies, or some form of obstruction to fluid 
flow. For example, in Hagen-Poiseuille pipe flow the fluid 
is confined to the interior of the pipe. Additionally, due 
to wall friction the fluid’s velocity at the wall is Zero. 
In a PHOENICS flow simulation problem, obstructions to 
fluid flow must also be introduced. Ina jet engine test 
facility model one must be concerned with obstructions to 
flow as a result of walls, floors, baffle material, turning 
vanes, and the jet engine itself. To account for the 


physical obstructions, PHOENICS makes use of the command 
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"CONPOR(R.R, TYPE,XF,XL,YF,YL,ZF,ZL)". The real number 
"R.R" normally has a value from 0.0 to 1.0. This value 
tells PHOENICS how much obstruction to flow is to be 
introduced by a surface or volume. For maximum obstruction 
(blocked flow), "R.R" would be 0.0. The "TYPE" argument 
specifies a surface or volume. For example a surface might 
be the inside of the test bay ceiling. A volume might be 
the geometric body of the jet engine. For a surface 
obstruction the "TYPE" Arguments "WEST, EAST, SOUTH, NORTH, 
LOW, OR HIGH" (the faces of an individual cell) are used. 
For a volume obstruction the type "CELL" is used. The 
arguments "XF, XL, YF, YL, ZF, and ZL" identify the first 
and last cells of the obstruction domain in the X, Y, and Z 
directions respectively. For example "CONPOR(0.0, HIGH, 
1,3,1,1,3,3)" tells PHOENICS to introduce complete 
obstruction to flow through the "HIGH" (Z direction) 
surfaces of the group of cells starting at IX=1, IY=1, and 
IZ=3, and ending at cells IX=3, IY=1, and IZ=3 
respectively. Each cell has an IX, IY, and IZ identity. 
Figure C.5 illustrates "HIGH" surface obstruction to flow. 
The cross-hatched surface blocks Z direction flow. "CONPOR 
(0.0,CELL,1,3,1,1,1,3)" would introduce maximum obstruction 
to flow through the surfaces of the volume defined by these 
coordinates. In Figure C.5, the "CELL" volume outlined 
with bold lines prevents fluid flow perpendicular to any of 


mes Surfaces. 
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PHOENICS also has the ability to introduce wall or 
surface "FRICTION" which brings about zero magnitude 
velocity components on a cell’s surface and calculates the 
friction coefficient. For the above "HIGH" example, the 
magnitude of flow in the Z direction was set to zero on the 
surfaces of several cells. To also set the X and Y 
velocity components to zero on the same "HIGH" surfaces, 
the CONPOR command requires a slight modification. 
Utilizing a "NEGATIVE" argument with the CONPOR command 
accomplishes the desired action. For example 
"CONPOR(0.0,HIGH,1,3,1,1,3,-3)" causes the X, Y, and Z 
velocity components on the "HIGH" surfaces of that 
particular group of cells to be zeroed. Through the use of 
the negative argument convention, velocity components on 
the external surfaces of a volume (cell) can be zeroed, and 


friction introduced into the model. 


F. SOURCE AND BOUNDARY CONDITIONS 

In any PHOENICS flow simulation problem description, 
the user’s Ql program must introduce source and boundary 
conditions to the EARTH program. Through the use of 
PHOENICS keywords, source and boundary conditions may be 
identified as energy, momentum, or mass flow. Energy 
sources are identified with the enthalpy variable "H", 
momentum by the velocity variables "U, V, or W", and mass 


flow utilizes the pressure variable "P". By using the 
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notation "H1" a fluid’s first phase enthalpy would be 
identified while "H2" would refer to second phase enthalpy 
for a two phase flow simulation. Information about inflow 
and outflow of any variable in the system domain is 
conveyed to the EARTH Program by way of the PHOENICS 
combination "PATCH-COVAL" statements. The PATCH command is 
used to specify the location at which a source or boundary 
condition is to be applied and to also specify the "TYPE" 
factor or what the PATCH command is to do to the system 
domain. The COVAL command is used to specify the 
coefficient and value which combine to determine the 
magnitude of the dependent variable over a PATCH. For 
example, "PATCH(JETOUT,LOW,1,1,5,6,8,8,1,1); 

COVAL(JETOUT, P1,FIXFLU,100.0); 
COVAL(JETOUT,W1,ONLYMS,500.0) ;COVAL(JETOUT,H1,ONLYMS, 
200000.0)" specifies a 100 kg/sec, steady-state, fixed-flux 
mass source. The mass source, named JETOUT, begins at the 
low faces of the cells IX=1, IY=5, and IZ=8 and ends at 
IX=1, IY=6, and IZ=8. The mass source has a Z direction 
velocity of 500 m/sec and an enthalpy of 200,000 J/kg. The 
COVAL coefficient "ONLYMS" specifies that only convection 
flow effects are influential in the transport of enthalpy 
from the cell faces to the cell centers. That is, 
diffusive effects are zero. Appendix 2 of the PHOENICS 


User’s Manual (Ref. 5) lists the various "TYPES", 
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"COEFFICIENTS", and "VALUES" that may be used with the 
combination "PATCH-COVAL" statements. 
G. INITIAL CONDITIONS 

In order to obtain the solution of a steady state 
problem, PHOENICS starts the iterative solution with an 
initial field of the dependent variables. The default 
initial fields for all dependent variables is zero over the 
entire physical domain. However, convergence to the true 
solution can be speeded up by the user specifying initial 
field conditions via "FIINIT" (Field Initialization) 
statements over the entire operational domain or 
"PATCH-INIT" statements over subdomains of the system. For 
example "FIINIT(H1)=AMOUNT" initialize a first phase 
enthalpy field to a user-defined value "AMOUNT". 
"PATCH (GASIN, INIVAL, XF, XL, YF, YL, ZF, ZL, TF, TL) ; INIT(GASIN, Pl, 
COEFFICIENT, VALUE)" initializes a first phase pressure 
field for the user-defined PATCH "GASIN". xXF,XL,YF,YL,ZF, 
and ZL define the geometric region of the PATCH called 


"GASIN". 


H. MAXIMUM/MINIMUM VALUES 

At the beginning of a flow simulation problem iterative 
process, a poorly defined dependent variable field value 
may have an adverse affect on problem convergence. At 
times, dependent variables "STRAY" outside of a physically 


meaningful range and cause numerical divergence. To help 
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minimize divergence of the solution PHOENICS utilizes the 
"VARMIN" and VARMAX" commands. These commands are used to 
set upper and lower bounds on any dependent variables. 
"VARMIN" (VI) =SLOW; VARMAX (W1)=FAST" sets a lower limit of 
"SLOW" to V1 and an upper limit of "FAST" to W1. 
I. UNDER-RELAXATION FACTORS 

A feature of PHOENICS which greatly assists problem 
convergence is the application of "UNDER-RELAXATION" 
factors to flow simulation dependent variables. This 
practice is most useful on an operational domain comprised 
of ill-proportioned cells (very large cells adjacent to 
very small cells). During program execution, EARTH’S 
dependent variable correction equations may develop 
enormously large coefficients between adjacent cells. The 
coefficient disparity often contributes to solution 
divergence. With selective "UNDER-RELAXATION" applied to a 
flow simulation program, PHOENICS introducers "FALSE-TIME- 
STEP" values to select dependant variables thereby 
preventing large dependent variable changes in magnitude 
per sweep. The "FALSE-TIME-STEP" brings equation 
coefficients closer together "LESS ABRUPTLY" which promotes 
problem convergence. "UNDER-RELAXATION" for pressure, 
velocity, and enthalpy components is frequently required 


for steady state problems. 


Za 


Je PROGRAM RUN PROCEDURES 

The Denmark test facility model was analyzed on the 
Naval Postgraduate School IBM-3033 mainframe computer 
system. Program run procedures involve interactive 
terminal sessions and a batch submission. The necessary 
flow simulation Q1 program and Job Control Language (JCL) 
programs reside in the User’s mainframe storage area. All 
flow simulation programs are submitted for execution by 
means of "EXEC" (execute) and JCL files. The EXEC files 
are activated by the keywords "COPYQ1", "RUNSAT", and 
"RUNEAR". 

The COPYQ1 EXEC file queries the user as to the name 
and file type of the flow simulation Q1 program source 
code. Afterwards the COPYQ1 EXEC file stores the source 
code and appropriate JCL in a temporary file. At this 
point, the user must ensure the above transfer of data was 
successful. This is accomplished through the user’s 
"VIRTUAL READER" file. An inspection of this file will 
reveal various condition codes which indicate the status of 
Gata transfer. Figure C.6 is an example of COPYQ1 EXEC and 
JCL files. 

The "RUNSAT" EXEC file extracts from the user’s 
mainframe storage area JCL which enables a Q1 program to 
EARTH program communications link. This JCL is also stored 
in a temporary file. Again a virtual reader file 


inspection must be made to ensure that a successful data 
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transfer was achieved. Figures C.7 and C.8 are examples of 
RUNSAT EXEC and JCL files. 

Finally after successful COPYQ1 and RUNSAT data 
transfers, the flow Simulation problem is ready for 
execution. The "RUNEAR" EXEC file is called, Necessary 
mainframe storage is allocated and the PHOENICS system is 
activated. Figure C.9 is an example of RUNEAR EXEC and JCL 
files. Upon program completion, PHOENICS output is sent to 
the user’s virtual reader file. After the PHOENICS output 
is analyzed, the user may discard the information or 
utilize a PHOENICS restart which allows the problem to 
continue from where it last ended. Upon completion of a 
PHOENICS run, output is stored on disk file "DFO9". To 
access this file and to continue the problem the PHOENICS 
code "RESTRT(A,B,....,Z)" must be utilized in the Ql 
program. "A,B,...,Z" represent the dependent variables 
which are to continue from their last value when the 
program ended. If the user desires all dependent variables 
to restart from their last values then "RESTRT(ALL)" can be 


used in the Q1 program. 


23 


IV. MODELING THE DENMARK TEST FACILITY 


A. MODEL DESCRIPTION 

The Denmark test facility model utilizes a Cartesian 
grid arrangement. Figure C.10 is a schematic diagram of 
the test facility in the Y-Z pane (side view). Figure C.11 
is the X-Z plane (top view) representation. The PHOENICS 
code allows for a plane of symmetry. For the Denmark test 
facility model the vertical Y-Z plane is one of symmetry. 
Therefore, in Figure C.11 only half of the Denmark test 
facility is represented, In this figure the jet engine is 
Slightly off-of-center from the Z axis for clarity only. 

To facilitate the Denmark test facility model geometry 
specifications, the system domain was divided into user 
defined zones. Each zone comprises a particular section of 
the Denmark test facility model. For example in Figure 
C.12, "YLBL1" is a user defined variable which specified 
the length of the first zone in the Y direction. This 
particular zone represents the length in the Y direction 
from the ground or floor of the Denmark test facility up 
2.0 meters. "ZLEXT2" represents the length in the Z 
direction of the atmosphere behind the test facility, and 
so on. Figure C.12 is a schematic diagram of the system 
domain zones which make up the vertical Y-Z plane while 


C.13 illustrates the zones in the horizontal X-Z plane. 
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The zone concept has proven to be very flexible. The 
number of cells in each zone may be specified by a user 
defined variable. For example, "XNCBL1" specifies the 
number of cells in the X direction for the length of 
operational domain "XLBL1". "ZNCEX2" specifies the number 
of cells in the Z direction for the length of operational 
domain "ZLEXT2" and so on. Grid refinement for a 
particular zone can be accomplished without any difficulty 
by changing one variable in the Q1 code. Figure C.14 is 
the method-of-pairs coding sequence which specifies the 
operational domain zones and numbers of cells used in the 
Denmark test facility model. Figure C.15 is an example of 
the EARTH output for the Q1 code in Appendix A. 

Program run time, program cost and EARTH output are 
directly affected by the number of cells in the model of 
the physical domain. By increasing the number of cells in 
the X, Y, or Z direction zones, the EARTH program requires 
more "SLAB-WISE" Computations thereby increasing computer 
CPU requirements and costs. By increasing the number of 
cells in the X or Y directions (i.e., NX or NY), the number 
of cells and therefore the number of calculations is 
increased. For example for a fixed NY, SAY 16, an increase 
of 1 cell in the X direction will result in 7 (number of 
dependent variables) times 16 = 102 additional equations 


(and unknowns). An increase of an additional cell in the Z 
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direction will result in 7 times NX times NY additional 
equations. For example for NX=9 and NY=16, increasing NZ 
from 28 to 29 will result in an additional 1,008 equations. 
Therefore, deciding to use a "COURSE" or "FINE" model of 
the physical domain is determined by the user’s time and 
financial constraints and the desired information from the 
flow simulation program. The Denmark test facility model 
has NX=9, NY=16 and NZ=28 for a grand total of 4,032 cells. 
Since each cell has 7 dependent variables, 28,224 linear 
algebraic equations must be solved iteratively to achieve 
problem convergence. This is a big job even for an IBM 
mainframe computer. 

For the purpose of this research, the Denmark test 
facility was modeled with a "COARSE" grid with the intent 
of achieving problem convergence prior to using a "fine" 
grid. Due to time constaints only one grid was used in 
this research. Therefore, problem convergence could not be 
identified as "GRID INDEPENDENT". 

Although the Denmark test facility model grid is 
considered "COARSE", many of the cells vary in size. The 
number and size of cells for a particular region of the 
grid was chosen based upon the expected activity ina 
region. For example, at the primary air inlet, one would 
expect primary cooling air to be entering at a relatively 
low velocity and in a laminar and uniform flow. Then as 


the primary air approaches the jet engine intake, fluid 
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velocities should be increasing in magnitude and changing 
in direction. At the outlet of the jet engine there exists 
a very hot and a very fast fluid flow mixing with entrained 
primary cooling air. In this particular region with a lot 
or activity, a less "COARSE" grid is more desirable in 
order to incorporate the mixing process into the model. As 
the exhaust gas enters the augmentor tube it mixes with the 
cooler entrained secondary air. Again, a less "COARSE" 
grid is necessary to better analyze the activity involved. 

Figure C.16 is a schematic diagram of the system domain 
grid for the X-Z plane while Figure C.17 represents the Y-Z 
plane. The two figures illustrate a less "COARSE" grid 
arrangement at the jet engine exhaust, augmentor tube 
inlet, and secondary air inlet. Again, grid refinement for 
any physical domain zone is accomplished by changing one 
user defined variable in the Q1 code. 

Although the Denmark test facility was modeled as a 
three-dimensional structure. The body fitted coordinate 
option of PHOENICS was not utilized. As such, the 
cylindrical jet engine, the augmentor tube, and the exhaust 
stack turning vanes were modeled as rectangular cubes. The 
jet engine was specified as a solid block with a hot mass 
source at one end and a mass sink at the other end. The 
jet engine was modeled with CONPOR and PATCH-COVAL 


statements. 
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The Denmark test facility walls, floor, and ceilings 
were modeled utilizing CONPOR and PATCH-COVAL statements. 
The CONPOR statements were given a zero porosity value 
thereby preventing mass transfer across any solid 
structure. PATCH-COVAL statements were necessary for the 
walls, floor, and ceilings in order to introduce 
appropriate friction and set fluid velocity components to 


zero magnitude at the structure surfaces. 


B. MODEL BOUNDARY CONDITIONS AND TURBULENCE MODELS 

Boundary conditions for the Denmark test facility 
openings to the atmosphere were specified with PATCH-COVAL 
statements. Ambient pressure was set to 101,325 Pascals 
while ambient temperature was set to 295 K. 

To incorporate turbulence into the flow simulation 
problem, this research utilized the PHOENICS K-E turbulence 
model. The K-E model introduces the KE (turbulence kinetic 
energy) and EP (kinetic energy dissipation rate) dependent 
Variables into the problem. Turbulence modeling is not an 
exact science and the choice of the K-E model for this 
research was arbitrary. The PHOENICS user’s manual 
(Ref. 5) explains in detail the K-E turbulence model and 
three other turbulence models available to the PHOENICS 
user. 

Readers may refer to Appendix A, the annotated Q1 code, 


which explains in detail the user defined variables for the 
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Denmark test facility model. Columns 1 and 2 are reserved 
for the beginning of executable statements, while non- 
executable statements (comments) begin in column 3 or 


beyond. 
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V. ANALYSIS OF RESULTS 


A. OBTAINING A FINAL SOLUTION 

The original purpose of this research was to computer 
model the Royal Danish Air Force Gas Turbine Test Facility 
using the PHOENICS code. The computer model would then be 
used to compare field data obtained from the test facility 
with computer generated pressures, temperatures, and 
velocities. If the computer generated output compared well 
with the field data, one might suggest that the computer 
model was successful. However, if field data from a test 
facility was non-existent or perhaps very limited how then 
could the success of a computer model be measured? 

The field data from the Denmark test facility was 
limited in scope. Besides the mass flow described in 
Chapter II, the only other field data available are the 
test bay pressure and the velocity and temperature of the 
engine exhaust at one point in the augmentor tube and at 
one point in the exhaust stack. Since the field data from 
the Denmark test facility was limited, a comprehensive 
comparison of field data with numerical results was not 
possible. To investigate the EARTH output (numerical 
results) and conclude that the numerical results were 
representative of the test facility, this study examined 


the EARTH output from three viewpoints: common sense 
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approach, conservation equation balances, and problem 


convergence. 


B. COMMON SENSE APPROACH 
The Denmark test facility was modeled with five major 
regions or zones. Within each zone a common sense approach 
to what one would expect the gas pressure, temperature, and 
velocity conditions to be was taken. The five major model 
zones are: primary air inlet, test bay, secondary air 
inlet, augmentor tube, and exhaust stack. The EARTH output 
presented in Appendix B is the numerical results of only 5 
slabs in the major model zones. 
1. Primary Air Inlet 

The primary air inlet is comprised of the following 
cells: IX=1 to IX=8, IY=1 to IY=12, and IZ=2. Cool ambient 
air (295 K) is drawn into the test facility through the 
primary air inlet (test facility doors). Pressure drops at 
the test facility doors were very small with a maximum drop 
of about 12 Pascals. Air flow through the primary air 
inlet was predominantly in the Z-direction (W). Velocities 
ranged from 1 to 2 m/sec with an average of about 1.7 
m/sec. The horizontal X-direction velocity (U) and the 
vertical Y-direction velocity (V) components were 
negligible. 

With an average gas density of about 1.2 kg/m? at 


the test facility doors, a primary air mass flow rate of 
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about 135 kg/sec was calculated. Recall from chapter IV 
that the Ql program modeled only one half of the test 
facility due to the Y-Z plane of symmetry. Therefore, the 
135 kg/sec primary air mass flow rate is for one half of 
the structure. The calculated primary air mass flow rate 
for the entire structure would be about 270 kg/sec. Of 
this amount, about 170 kg/sec would be primary cooling air 
and the remaining 100 kg/sec flows through the engine. 
From this point on, references to mass flow rates will be 
for the entire test facility. The temperature at the test 
facility doors remained at about 295 K. At the test 
facility doors, one would not expect much activity and a 
laminar flow. The KE and EP turbulence dependant variables 
were negligible which indicated laminar flow. 
2. Test Bay 

The test bay comprises cells IX=1 to IX=8, IY=1 to 
TY=12, and IZ=2 to IZ=9. The jet engine which was modeled 
aS a 2000 K, 100 kg/sec sink~source is located at cells 
IX=1, IY=5 to IY=6, and IZ=6 to IZ=7. At 2000 K, the 
density of the engine exhaust gasS was 0.18 kg/m>. The jet 
engine provides the excitation which draws ambient air in 
through the primary air inlet, entrains primary cooling air 
into the augmentor tube from the test bay and entrains 
secondary cooling air into the augmentor tube from the 
secondary air inlets. As primary air moved through the 


test bay from the primary air inlet to the sink side of the 
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jet engine, its pressure and velocity changed. At the two 
cells (IX=1, IY=5 to IY=6, and IZ=5) before the engine 
inlet, about 1.5 meters upstream of the inlet, primary air 
pressure dropped about 40 Pascals. The primary air W- 
velocity component was about 3 m/sec at IZ=4, about 3 
meters upstream of the engine inlet. The horizontal U- 
velocity component was about -4.5 m/sec which indicated 
that air was being drawn into the engine from the side 
walls of the test bay. The vertical V-velocity component 
was about -4 m/sec above the engine and about 4 m/sec below 
the engine. This indicated that air was being drawn down 
from the test bay ceiling into the engine and at the same 
time being drawn up from the test bay floor. The turbulent 
parameters KE and EP were very small indicating laminar 
behavior at the engine inlet. The temperature at the 
engine inlet remained at about 295 K. 

At the discharge side of the engine (cells IX=1, 
TY=5 to IY=6, and IZ=8), significant activity was observed. 
The pressure increased about 2100 Pascals. In adjacent 
cells above/below and longside the engine at IZ=8, the 
pressure dropped about 410 Pascals. The pressure 
difference among the cells caused an expanding or 
diverging effect if the jet engine exhaust gas. The U- 
velocity component of the gas was about 30 m/sec while the 
V-velocity component was about 50 m/sec. The jet engine 


exhaust gas expanded out towards the rear wall of the test 
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facility as opposed to going directly into the augmentor 
tube. However, due to the large magnitude of the W- 
velocity component, about 430 m/sec, the resultant 
direction of the engine exhaust was primarily into the 
augmentor tube. Hot, fast moving gas at the engine outlet 
mixed with the cool, slow moving entrained primary cooling 
air. Due to this mixing, small turbulence resulted. this 
premise 1S supported by the magnitude of the KE and EP 
parameters in the region. KE had a maximum value of about 
110 J/kg, while EP had a maximum value of about 5400 
J/kg/sec. Due to this mixing of engine exhaust and primary 
cooling air, over a distance of about 1.5 meters, the 
exhaust gas was cooled from 2000 K to about 1830 K prior to 
entering the augmentor tube (IZ=10). Over the same 1.5 
meter distance, the exhaust gas W-velocity component 
decreased from about 430 m/sec to about 400 m/sec. 

The investigation of the test bay was localized to 
a small region around the jet engine. For the most part 
the activity in the rest of the test bay was uneventful 
(1.e. no turbulence and insignificant pressure, 
temperature, or velocity changes). 

The above observations in the test bay and at the 
inlet and outlet ends of the jet engine are realistic. One 
would expect primary air velocity components and pressure 
to change somewhat prior to going into the engine and more 


Significantly at the engine outlet. Additionaliy, due to 
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the primary air being drawn into the engine a low pressure 
area at the inlet makes good sense. While at the outlet 
side of the engine a high pressure is expected. At the 
engine outlet one would expect that the hot, fast engine 
exhaust colliding and mixing with the cool, slow entrained 
primary cooling air would result in some turbulence. 

3. Secondary Air Inlet 

As was mentioned earlier in chapter II, the 
secondary cooling air plays a significant role in the 
operation of a jet engine test facility. The secondary 
cooling air serves the purpose of providing additional cool 
(ambient) air for mixing with the hot engine exhaust. For 
the Denmark test facility the largest mass flow rate occurs 
through the secondary air inlets. The test facility data 
indicated a mass flow rate of about 340 kg/sec. In the 
structure model, secondary cooling air enters the augmentor 
tube at cells IX=5, IY=1 to IY=10, and IZ=10 to IZ=12. The 
secondary cooling air had an average U-velocity component 
of about -14 m/sec, which indicated that secondary cooling 
air was being drawn into the augmentor tube from outside 
the structure. 

A maximum pressure drop of about 480 Pascals at the 
secondary air inlet was in consonance with the relatively 
high U-velocity component. With an average air density of 
aeout 1.0 kg/m? at the secondary air inlet, a mass flow 


rate of about 420 kg/sec was calculated. Although the 
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secondary air inlet had a significant pressure drop and a 
relatively large U-velocity component (about -14 m/sec) it 
did not demonstrate turbulent behavior based upon the 
relatively low KE and EP parameters. The temperature at 
the secondary air inlet indicated a temperature rise of 
about 15 K. The results showed that warm air from the 
exhaust stack was flowing into the secondary air inlet 
region. 

With a primary cooling air mass flow rate of about 
170 kg/sec and a secondary cooling air mass flow rate of 
about 420 kg/sec, an Augmentation Ratio of about 5.9:1 was 
calculated. 

4. Augmentor Tube 

Three regions of the augmentor tube were chosen for 
investigation: The region IZ=10 to IZ=16, the slab IZ=20, 
and the slab IZ=22. At the beginning of the augmentor 
tube, IZ=10 to IZ=16, a large mass of primary and secondary 
cooling air collided and mixed with the hot, fast exhaust 
gas. As a result, Significant turbulence occurred. A 
maximum KE value of about 7,100 J/kg was observed while 
maximum EP was about 500,000 J/kg/sec. In the 8.5 meters 
that the engine exhaust traveled, it cooled from 2000 K to 
about 920 K. The engine exhaust W-velocity component 
decreased from about 430 m/sec to about 125 m/sec. The U 
and V-velocity components were considered negligible. 


Figures C.18 and C.19 are PHOENICS enthalpy (H1) and 
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velocity (W1) profiles for the cells IX=1, IY=1 to IY=10, 
and IZ=16. The reader is reminded at this point that an 
enthalpy value divided by a specific heat value of about 
1000 J/(kg K) will give a good temperature approximation 
within a cell. The profiles indicated minimum temperature 
and W-velocity values at the bottom (IY=1) and top (IY=10) 
of the augmentor tube and maximum temperature and W- 
velocity values at the centerline of the augmentor tube. 
These temperature and W-velocity distributions make good 
sense. Since the location IY=5 is at the height of the 
engine outlet, one would expect the engine exhaust to be 
hotter and to be moving faster here than at the top/bottom 
of the augmentor tube. Figures C.20 and C.21 are 
temperature and W-velocity profiles for the cells at IxX=1 
to IX=5, IY=5, and IZ=16. Although these profiles are at 
the same slab location (1Z=16) as figures C.18 and C.19, 
the examination of the variables is from the center of the 
augmentor tube to the side wall. These distributions 
indicated maximum values of the variables occurred at the 
center of the augmentor tube and minimum values at the side 
walls. These profiles also make good sense. In the cells 
at the augmentor wall (1IX=5, IY=1 to IY=10, and IZ=16), the 
secondary cooling air entering the augmentor tube has only 
mixed with the engine exhaust for about a 3 meter distance 


down the augmentor tube. 
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At augmentor tube location IZ=20 (about 22 meters 
downstream of the engine outlet), the total mass flow was 
examined. The W-velocity component ranged from about 25 
m/sec to 75 m/sec with an average of about 45 m/sec. The U 
and V-velocity components were considered to be negligible. 
With an average gas density of about 0.6 kg/m, the total 
mass flow rate was about 700 kg/sec. At augmentor tube 
location IZ=20, there was still some turbulence, which 
indicated continued mixing of the primary/secondary cooling 
air and the engine exhaust. The maximum KE value was about 
740 J/kg while EP maximum was about 10,400 J/kg/sec. The 
temperature ranged from about 520 K to about 635 K with an 
average of about 550 K. Figures C.22 and C.23 are 
temperature and W-velocity profiles for the cells at IX=1, 
ITY=1 to IY=10, and IZ=20. Figures C.24 and C.25 are 
temperature and W-velocity profiles for the cells at IX=1 
to IX=5, IY=5, and IZ=20. A more uniform distribution of 
temperature and W-velocity was observed across the slab at 
IZ=20 than at slab IZ=16. This is reasonable since the 
total mass of primary and secondary cooling air and engine 
exhaust has mixed for a distance of about 22 meters 
downstream of the engine outlet. 

In the augmentor tube/exhaust stack transition 
region (IZ=22), the test facility field data indicated the 
W-velocity component on a centerline with the jet engine 


was about 60 m/sec and the average temperature was about 
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650 to 675 K. The numerical results indicated a W-velocity 
component of about 63 m/sec and a temperature of about 

615 K. Ina distance of about 4 meters from IZ=20 to 
IZ=22, turbulence decreased by about 50 percent. KE was 
about 400 J/kg while EP was about 4,500 J/kg/sec. 

5. Exhaust Stack 
Field data from the Denmark test facility indicated 

a temperature rise of about 20 K from the bottom of the 
exhaust stack to the area under the exhaust stack baffle 
material. In the numerical model a temperature drop of 
about 50 K was observed. The Denmark test facility was 
modeled without baffle material at the top of the stack. 
Gas velocities at the top of the computer model stack were 
for the most part in the V-direction. The V-velocity 
ranged from about 10 to 45 m/sec with an average of about 
20 m/sec. Turbulence at the top of the exhaust stack was 
considered to be negligible. Average temperature at the 
top of the exhaust stack was about 550 K. Figures C.26, 
C.27, C.28, and C.29 are temperature and V-velocity 
profiles for the top (IY=13) of the exhaust stack. Figures 
C.26 and C.27 are for the cells IX=1, IY=13, and IZ=21 to 
IZ=26 (along the Z axis). Figure C.26 indicated highest 
temperature values towards the rear (IZ=25) of the exhaust 
stack. Additionally, Figure C.27 indicated that the 
V-velocity fluctuated (low to high to low ...) from IZ=21 


to IZ=26. It is suggested that the manner in which the 


CA 
(0 


turning vanes were computer modeled caused a lack of 
uniformity in the temperature and V-velocity distribution 
for this stack region. 

Figures C.28 and C.29 are the temperature and V- 
velocity distributions for the cells at IX=1 to IX=5, 
TY=13, and IZ=23 (from the center of the exhaust stack out 
to the side wall). These figures indicated fairly uniform 
temperature and V-velocity profiles throughout the cells in 


this region. 


C. COMPUTER DATA AND FIELD DATA COMPARISON 

The difference between field data results and numerical 
results can be attributed to several things: baffle 
material (not included in the numerical model), turning 
vanes (crudely modeled as rectangular cubes), turbulence 
model (modeling is not an exact science), and model grid 
(grid independence not verified). 

The computer model total mass flow rate of about 700 
kg/sec is about 15 percent larger than the 600 kg/sec total 
mass flow rate of the Denmark test facility. It is 
suggested that not including the baffle material, in the 
computer model primary and secondary air inlets and the 
exhaust stack, resulted in larger numerical mass flow rates 
than the test facility mass flow rates. Test facility 


baffle information was not available for this research. 


Field data indicated a maximum pressure drop in the 
test bay of about 500 Pascals. The computer model 
calculated a maximum pressure drop of about 410 Pascals in 
the cells adjacent to the jet engine outlet. 

AS was mentioned earlier, the exhaust gas temperature 
and W-velocity component in the augmentor tube/exhaust 
stack transition region compared very well with the field 


data. 


D. CONSERVATION EQUATION BALANCES 

The PHOENICS output code enables the investigator to 
determine the degree conservation of mass has been 
achieved. The PHOENICS output symbol for a mass flow 
source is "R1". By an algebraic summing of the mass flow 
sources, a conservation of mass flow calculation can be 
made. For a total accounting of all the mass flow sources 
in a computer model, the algebraic sum of the sources 
should be zero. For the Denmark test facility this sum was 
equal to 0.02 kg/sec. If this amount is compared to the 
100 kg/sec mass flow through the jet engine there is a 0.02 
percent error. Compared to the 700 Kg/sec mass flow 
calculated for the augmentor tube the percent of error is 
reduced to 0.003. Based on this analysis, a mass balance 


Or a continuity of flow exists in the Denmark test facility 
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model. A mass balance or a continuity of flow is a 
necessary but not sufficient condition for problem 


convergence. 


E. CONVERGENCE 

Chapter III described the iterative process which 
adjusts the values of the dependent variables during the 
solution process. During the process, imbalances occurred 
between the left and right hand side of the dependent 
variable discretized algebraic equations. These imbalances 
are called residuals. Problem convergence would be 
obtained when all residuals were zero (i.e. the dependent 
variable discretized algebraic equations were in balance). 
The residual goes to zero as the number of sweeps goes to 
infinity. With numerical solutions, users must employ a 
variety of methods to determine problem convergence since 
residuals can not decrease to zero, as a practical matter. 

With the PHOENICS output (Appendix B), this research 
utilized a local and a global analysis of the dependent 
variables to determine problem convergence. A local or 
"spot value" analysis can be accomplished with the PHOENICS 
output by designating cells to be monitored. Since the 
test facility field data referenced the exhaust gas 
temperature and velocity for a location in the augmentor 
tube/exhaust stack transition region, the computer model 


spot checked the dependent variable values at the same 
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ied oa 


location (the cell at IX=1, IY=5, and IX=22). From the 
initial program run to the final sweep (about 10,000), the 
spot values were closely monitored and compared. Local 
convergence at this cell occurred when little or no change 
in the dependent variable values was observed from sweep to 
sweep. Local or spot value convergence of dependent 
variables for one cell is a necessary but not sufficient 
condition for problem convergence. For this reason, the 
results at "ALL" cells were spot checked during the 
solution process. 

During the problem solution process, the PHOENICS code 
sums the dependent variable whole-field residuals. An 
examination of the whole-field sums provides a means of 
global investigation to determine problem convergence. If 
whole-field residuals decrease towards zero, with 
increasing iterations, problem convergence is suggested. 
However, if the residuals increase, with increasing 
iterations, problem divergence is occurring and some 
adjustment of the PHOENICS relaxation parameters is 
required. Global analysis is also accomplished by 
examination of the energy balance of the system. The 
dependent energy variable enthalpy (H1) residual is first 
Givided by 1,000,000. This result is then divided by the 
absolute value of the sum of the enthalpy (H1) net sources. 
If the quotient decreases to less than about 0.05 (5 
percent error), problem convergence is suggested. If the 
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quotient continues to increase, problem divergence is 
occuring and PHOENICS relaxation adjustment is again 


required. 


F. SUMMARY 

Although there was limited test facility field data to 
compare with numerical results, it is suggested that the 
computer modeling of the Denmark gas turbine test facility 
with PHOENICS code was successful. The field test data 
compared well with the numerical results. During the 
common sense examination, the conditions that existed in 
the five major regions of investigation made good sense and 
were indicative of what one might expect in a gas turbine 
test facility. The total mass flow continuity check 
indicated less than 0.003 percent error. Finally, local 
and global convergence was indicated. A system energy 


balance was within a 5 percent error. 
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VI. RECOMMENDATIONS 


A. INTRODUCTION 

It is believed the PHOENICS code is a useful tool for 
the solution of gas turbine flow simulation problems. 
However, since the Denmark test facility has not "FINELY" 
modeled the following computer model enhancements are 
suggested: BFC, PHOTON, alternate turbulence models, grid 


refinement, and baffle modeling. 


B. BFC 

Body fitted coordinates (BFC) are necessary to 
incorporate the cylindrical and curved surfaces of the jet 
engine, augmentor tube, and turning vanes. With the proper 
geometric shapes in the test facility model, mass flow in 
the augmentor tube inlet, the secondary air inlets, and the 
exhaust stack will be more representative of the actual 


test facility. 


C. PHOTON 

The implementation of the PHOENICS PHOTON option will 
greatly enhance numerical results investigation. A two or 
three dimensional representation of a dependent variable 
field is much more useful than a cell by cell examination 


of numbers. 
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D. ALTERNATE TURBULENCE MODELS 

One or more of the alternate PHOENICS turbulence models 
should be implemented and investigated. Since it was 
suggested that turbulence modeling is not an exact science, 
verifying the effect of the turbulence model on the 
solution is crucial. It may be that the turbulence model 


selected does not greatly effect the solution. 


E. GRID REFINEMENT 

Increasing the number of cells (grid refinement) in the 
system will enhance the problem solution process. During 
the slab-sweep iterative process, the dependent variable 
discretized algebraic equations will be solved ina "LESS 
ABRUPT" manner which will promote problem convergence. 
Grid refinement is also necessary to verify the solution is 


grid independent. 


F. BAFFLE MODELING 

Finally, by utilizing additional PHOENICS CONPOR 
commands, baffle material could be modeled into the 
primary/secondary air inlets and the exhaust stack. By 
experimenting with different values of porosity, the total 
mass flow rate of the computer model could more closely 
represent the total mass flow rate of the actual test 


facial tye 
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APPENDIX A 


Q1 CODE 
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TALK=F 3RUN( 1,1) 


* 
9 


* 
% 


HICH HII HEE FEE FEIE FEE HE HE FEFE FE FEE HE FE HE HEHE HEHE FE IEFE DE HE HE HEFEHE FEE HE FEES PEFEHE TEFE HE FEFE HEHE EFEFE FETE DEKE HE IEE 


% 
% 


* KK KK KK KK KK KK KK KK KK KK KK KK KK KK K KK KOK OK OK OK OK OK OK OK KOK OK OK OK ORK OX 


DENMARK JET-ENGINE TEST FACILITY COMPUTER MODEL 
BE IEE IE HH IE IEE AL IE HE AEE KSEE HIE AHI IERIE IE HE IIE EH EI HIE IH HIE HEHE 


XLBLI 

XLMIDI 
XLMID2 
XLMID3 
XLBL2 

XLBL3 

XLEXT1 
XLEXT2 
*NCBL1 
XNCMD1 
XNCMD2 
XNCMD3 
XNCBL2 
XNCBL3 
XNCEX1 
XNCEX2 
YeEnID 

YLBLI 

WEEE? 

YLEXT1 
YLE*h2 
YLBL3 

YEEATS 
YLBLS 

YNCMID 
YNCBLI 
YNCBL2 
YNCEX1 
YNCEX2 
YNCBL3 
YNCEX3 
YNCBLG 
CLE Xie 
ZUING? 
2Eihihe 
eeINts 
ZLINTS 
ZLINTS 
ZLINT6 
CLEATZ 
ZNCEX1 
ZNCIN1 
ZNCIN2 
ZNCIN3 
ZNCING 


NON-PHOENICS USER-DEFINED VARIABLES 


THE FOLLOWING VARIABLES ARE USED TO DESCRIBE THE 


PHYSICAL CHARACTERISTICS OF THE DENMARK TEST 
THE FIRST LETTER OF THE VARIABLE 
SPECIFIES THE DIRECTION IN THE COORDINATE 


FACILITY. 


SYSTEM. 


THE FOLLOWING CONVENTIONS ARE USED TO DESCRIBE 
A PARTICULAR ZONE OF THE MODEL: 
LENGTH IN A COORDINATE DIRECTION 


E 
MID 
BL 
EXT 
NC 
INT 


EXAMPLE: 


MIDDLE OF SAY THE 


AUGMENTOR TUBE 


BOUNDARY LAYER 

EXTERIOR ZONE OF STRUCTURE 
NUMBER OF CELLS IN A DIRECTION 
INTERNAL ZONE OF STRUCTURE 


XLBLI 


YNCEX2 


ZLINTS 


THE LENGTH OF THE FIRST 
ZONE IN THE X DIRECTION. 


THE NUMBER OF CELLS FOR THE 
SECOND EXTERNAL ZONE IN THE Y 
DIRECTION. 


THE LENGTH OF THE FIFTH 
INTERNAL ZONE IN THE 2 
DIRECTION. 


A48 


% 
% 


KK KK KKK KK OK OK OK KK KK OK OK KK KK OK OK OK OK OK KK KK KK KK KOK KK KK KK KK OK OK OK 


* ZNCINS 

* ZNCING 

* ZNCEX2 

* 

% 

* TMPJET TEMPERATURE OF JET EXHAUST (XK) 

* CPJET SPECIFIC HEAT OF JET EXHAUST (J/KG#K } 
* PRSJET PRESSURE OF JET EXHAUST (N/M¥2 ) 

* GASCON GAS CONSTANT OF JET EXHAUST (N¥M/KGK ) 
* MASJET MASS FLOW RATE OF JET EXHAUST (KG/SEC ) 
* VELJET VELOCITY OF JET EXHAUST (M/SEC) 

* RHOJET DENSITY OF JET EXHAUST (KG/M*3 } 

* 

* TMPAMB AMBIENT AIR TEMPERATURE (K)} 

* PRSAMB AMBIENT AIR PRESSURE (N/MX¥*2 } 

* CPAMB SPECIFIC HEAT AMBIENT AIR (J/KGK ) 

* 

* 

2% 

“ 


* KKK KK KK OK KOK KK OK K OK KK OK OK 


HHH HIE IIE ESE HIE FEE IE HE EIEIO EIEIO SEE SE IE HE IEEE ESE JE SESE SE SESE IESE DE JE SE SE HE SESE IE FE FE FE SEIEDE SEE ESE DE HE 


* 
* 
* 


* 
* 
* 


HH IIE H HEHE IE HEHE HEE HE EE IEE HE EIEIE SESE HE IEE IE IESE IESE SE HE ESE IEE IESE SE EEE HE HEE IESE ESE JE IESE HE SESE IEEE EIEN JE HE HEE HE HE 


* 


* VARIABLE DECLARATIONS 
* 


* 
* 
* 


* 

* PHOENICS ALLOWS 45 USER DEFINED REAL-NUMBER VARIABLES 

% 
REAL(XLBL1,XLMID1,XLMIOZ2,XLMIO3,XLBL2 ,XLBL3,XLEXT1,XLEXT 2} 
REAL(YLMID,YLBL1,YLBL2,YLEXT1,YLEXT2,YLBL3 ,YLEXT3 ,YLBL } 
REAL(ZLEXT1,Z2LINT1,CLINT2 ,ZLINT3 ,ZLINT4¢,ZLINTS»,ZLINT6 ,ZLEXT 2 3 
REALCTMPJET ,CPJET ,PRSJET ,GASCON ,MASJET , VELJET , RHOJET 3} 
REAL( TMPAMB ,PRSAMB ,»CPAMB ,DELT,FACT ) 
REAL( KEINJ,EPINJ,KEINA,EPINA ) 

x 

# PHOENICS ALLOWS 45 USER DEFINED INTEGER VARIABLES 

x 
INTEGER( XNCBL1,XNCMD1 ,XNCMD2 »XNCMD3 ,XNCBL2 ,XNCBL3 »XNCEX1 ,XNCEX2 ) 
INTESERCYNCMID,YNCBL1,YNCBL2 ,YNCEX1 , YNCEX2 » YNCBL3 »>YNCEX3 » YNCBL4 ) 
INTEGER( ZNCEX1 ,ZNCINI »ZNCIN2Z ,ZNCINS »ZNCING »ZNCINS ,ZNCING » ZNCEXZ 3} 

3% 

2% 

* 


* 
*% 
* 


* 


* 
* 
* 


* 


x INITIALIZATION OF USER-DEFINED VARIABLES 
* 


* 


Ag 


* 
* 
* 


* 


~ X DIRECTION LENGTHS (METERS) 
* 

XLBL1=0.5 

XLMID1=0.5 

XLMID2=0.5 

*XLMID3=0.5 

XLBL2=0.5 

XLBL3=5.4 

MUEXEI=2.0 
3e 
x NUMBER OF CELLS IN X DIRECTION 
% 

XNCBL1=1 

XNCMD1=1 

XNCMD2=1 

XNCMD3=1 

XNCBL2=1 

XNCBL3=3 

XNCEX1=1 
* 
x Y DIRECTION LENGTHS (METERS) 
% 

YLBL1=2.0 

YLMID=1.0 

YLBL2=2.0 

YLEXT1=3.85 

YLEXT2=2.15 

YLBL3=6.0 
* 
_ NUMBER OF CELLS IN THE Y DIRECTION 
* 

YNCBL1=4 

YNCMID=2 

YNCBL2=4 

YNCEX1=2 

YNCEX2=1 

YNCBL3=3 
* 
- Z DIRECTION LENGTHS (METERS} 
x 

ZLEXT1=2.0 

ZLINT1=23.5 

ZLINT2=3.2 

ZLINT3=1.5 

ZLINT4=6.0 

ZLINTS=16.3 

ZLINT6=8.4 

ZULENT2=40 0 
ad 
x NUMBER OF CELLS IN THE Z DIRECTION 
* 

ZNCEX1=1 

ZNCIN1=5 

ZNCIN2=1 

ZNCIN3=2 


90 


* 


* 


*K 


* 


* 


ZNCING=6 
ZNCINS=5 
ZNCIN6=6 
ZNCEX2=2 
# 
* 


* 
* 
* 


HHH HEH HEH HE IEE HICH FEE FETE FE IE IE FE DE FE FE FEDE DE HE HE FEEL IE DE HEHEHE FE IE HE IESE HE HE FE FE HE HE ETE HEHEHE HEHEHE HE HEHE HEHE HEE HE HE HEHEHE 
HE FEI HEHEHE HE HE IE FIO HE IEE HE HEHE FETE HE HEHE HIE HE HEHE HE IEEE HE HE HEHE FEE HE HE HEHEHE HE EE EE SH HEHE HEH HEE HE HEH HH HEHE 


* 


* GROUP 1. 
* 


RUN TITLE AND OTHER PRELIMINARIES 


* 
* 
* 


HE HIKE HH HEHEHE HIE HE IEE HE EE HEHE FEE HE HEH EE HEE IE HE HEHE EE EE DIESE HE EEE FE FE HEHEHE EH IEE HE HEI IE FEE HE HEHE IESE 


* 


TEXT( DENMARK JET-ENGINE TEST FACILITY MODEL) 


* 
* 
* 


* 


* 
* 
* 


HHH KH IIH He EE HE HEHE HE IEE HE EEE HE EIEIO SE HE FEE IE FE HE FE HE HE IESE HE HE FE HE ESE HE HE FEE FEE HEHE EE IEE HEHE FEE HE HE HEHE HE 
HEH HEH HAH HH HEHEHE HE HEH HEHEHE HEH HEHEHE HEHE HH IH HEHE HEHEHE HEHE HEHE HEHE HEHEHE HI IOI HE HEE HEHE HEHE HE HEHEHE HE 


GROUP 2. 


* KK K KK K 


TRANSCIENCE 3 TIME-STEP SPECIFICATION 


DEFAULT TAKEN: STEADY STATE 


* * KK KK K 


HEH HEF HEH HHH HEH HH IIE HEH HEH HE HEH HIE IE HE HEE HEHE HEHEHE HEHE HEHEHE EH HEHE HEHEHE HEHE IEE HEHE HEHE HEHE HE IIE HE HE MK HK 
HEHEHE HEH IETHER HE HAH HH HHH H HHH 


* 


* GROUP 3. 
* 


A-OIRECTION GRID SPECIFICATION 


* 
* 
* 


HE FEIE HE HEHEHE EEE HEHE HEHE HE HEHE HE HEE EE EH EE EE HE HEE HEHE EE HEHE HEE HE HEHE HEHE FEE HE HE HEHE IIE HE HE HEHEHE HII HHH HK 


* 


* 


* NA = THESNUMBER CReCELLS IN BE» DIRECTION CINTEGER ) * 
* MULAST = DISTANCE IN THE X DIRECTION (METERS) * 
* * 
NX=XNCBL1+XNCMD1+XNCMD2 +XNCMD3 +XNCBL2+XNCBL3+XNCEX1 
x» 
XULAST=XL3BL1+XLMID1+XLMION2+XLMIDO3S+XLBL2+XLBL3+XLEXT1 
* * 
* THE METHOD-OF-PAIRS HAS BEEN USED TO DEFINE THE NUMBER * 
* AND SIZE OF CELLS. THE XFRAC,YFRAC, AND ZFRAC SPECIFY THE X ,Y>* 
* Bin) 2 COORDINATES RESPECTIVELY. manne COD INDEX SPECIFIES THE * 
* NUMBER OF CELLS IN A DIRECTION AND THE EVEN INOcxX SPECIFIES THE* 
* SIZe OF EACH CELL. THE NEGATIVE GUANTITY FOR THE FIRST OOD * 
* INDEX IN EACH GROUP IS THE RESERVED WORD WHICH SPECIFIES * 
* "METHOD-OF -PAIRS" * 
* * 


XFRAC(1 )=-XNCBLI; 
XFRAC(3 )=XNCMD1; 
XFRAC(S J=XNCMD 23 
XFRAC( 7 )=XNCMD33 
XFRAC(9)=XNCBL2; 
XFRAC( 11 )=XNCBL3 5 


XFRAC( 2 J=XLBLI/( XNCBL1¥XULAST ) 
XFRAC(4 J=XLMIDI/( XNCMD1L¥XULAST ) 
XFRAC( 6 )=XLMID2/( XNCMD2*XULAST ) 
XFRAC(8 J=XLMID3/(XNCMD3*XULAST } 
XFRAC(10 )=XLBL2/( XNCBLZ*XULAST ) 
XFRAC(12 )=XLBL3/( XNCBL3*XULAST ) 


ol 


XFRAC (13 J=XNCEX1 $5 XFRAC (14 )=XLEXT1/( XNCEX1*XULAST ) 


# * 
* * 
* * 


HH IEE HE IEE HEHE EIEIO IE HEA IE HEHE EEE HE IEE EH EIEIO IE FE IEE FE FEE HE FEHE FE HEHE FEE HEE HE FE HE SEE HE EEE FEE EE IE ESE SEE HE IEE HE 
FFE IEE FEE IE FE FEE IE HE FE HE IE IEE EEE EE FEE IE FE IEE FE FEE EE HE FE FEE HEIL FE FE IEE HE HEI EE IEE HE HEHEHE EH IEEE FEIN 


# * 
* GROUP 4. Y-DIRECTION GRID SPECIFICATION * 
x 
BEIM HIE BEE IE IE HEHE HEME HOE HIE HEISE IH HEHE EHO HIME EI IEE ITE TI OL HEE TE IE 
* * 
* NY = THE NUMBER OF CELLS IN THE Y DIRECTION (INTEGER) * 
* YVLAST = DISTANCE IN THE Y DIRECTION (METERS) * 
- x 


NY =YNCMID+YNCBL1+YNCBL2+YNCEX1 +YNCEX2+YNCBL3 
YVLAST=YLMID+YLBL1+YLBELZ+YLEXTI+YLEXT2+YLBL3 


* * 
* METHOD OF PAIRS * 
x x 

YFRAC(1)=-YNCBL1} YFRAC(2 )=YLBL1/( YNCBL1*¥YVLAST } 

YFRAC(3 )=YNCMID; YFRAC(4 J=YLMID/( YNCMID*YVLAST ) 

YFRAC(5 J=YNCBL2; YFRAC(6 J=YLBL2/( YNCBL2*YVLAST ) 

YFRAC( 7 J=YNCEX13 YFRAC(8 J=YLEXT1/( YNCEX1*YVLAST ) 

YFRAC( 9 )=YNCEK235 YFRAC(10 )=YLEXT2/( YNCEX2*YVLAST ) 

YFRAC(1L1J=YNCBL33 YFRAC(12 J)=YLBL2/( YNCBL3*YVLAST ) 
* * 
ie x 
x * 


HHH IIH HAH HII HII IE HEHEHE HEHE IEE EE HEHE HEHE HEHE IH HEHE HEH HEHEHE HE HEHE EH IE HIE IE HEHE EH HEHE IIE HHH HE HEHE HE 
HH HHH HAHAH AH HAH HK AHHH HHH IH HHH HII HEHE HEH HEH IH HEA IH HH HEHE HHI HEI HHH HHH AHH RAKE 


* x 
* GROUP 5. Z-DIRECTION GRID SPECIFICATION * 
* * 
HHH HHH HEE HEHE HEHEHE HEHEHE HEHEHE HEHEHE HEHEHE HE IEEE IE IESE HEHE HEHE HIE HEHEHE HE HEHE HEHEHE HE HEHE 
x 
* NZ = THE NUMBER OF CELLS IN THE 2 DIRECTION (INTEGER ) * 
* ZNLAST = DISTANCE IN THE 2 DIRECTION (METERS) * 
~ a 

NZ=2NCEX1+ZNC IN1+#ZNCIN2Z+#ZNC IN3 +ZNCING+ZNCINS+ZNCING +ZNCEX2 

CRLAST=CLEXTL4+ZLINT1+ZLINT2+ZLINT3+2ZLINT4¢+ZLINTS+ZLINT6+ZLEXT 2 

ZFRAC(1)=-ZNCEX1; ZFRAC( 2 )=CLEXT1/( ZNCEXL*ZWLAST ) 

ZFRAC(3 J=ZNCINI 5 ZFRAC(4 J=ZLINTI/( ZNCIN1*2WLAST ) 

CFRAC(5 )=2NCINZ, ZFRAC( 6 JSZLINT2/( ZNCINZ*ZWLAST ) 

CFRAC( 7 J=Z2NCING 5 ZFRAC( 8 J=ZLINTS/( CNCIN3Z*ZWLAST ) 

ZFRAC( 9 )J=ZNCING; ZFRAC( 10 J=ZLINT4/( ZNCING*ZWLAST ) 

ZFRAC( 11 J=Z2NCINS; ZFRAC( 12 )=ZLINT5/( CNCINS*ZWLAST ) 

ZFRAC(13 JEZNCING 3 ZFRAC(C14 JEZLINTO/( ZNCING*ZWLAST ) 

ZFRAC(15 J=Z2NCEX23 ZFRAC( 16 J=ZLEXT2/( 2NCEX2*2ZWLAST ) 
> x 
* x 
x x 


HH HHI IE HEHE HIE FEE IE HE FEIE HE FE IE FH EIE HE HEHEHE HEHEHE EEE HEHE PEE HEHEHE HEE HE HEE EE EEE EOE FE HE EIEN HEE HEHEHE HEHE HE HEHE 
HHH HAF HHH IIH HEHE IE HE HEHEHE HE HIE HEHEHE HII HEHE IEE HIE HEHE IEE HE HEHEHE HE IER HEH IE HEHE HEHE HEHE FEE HE HEHE 
* * 


* GROUP 6. BODY-FITTED COORDINATES OR GRID DIRECTION * 


* (NOT IN USE) * 
FE FEI IEE IEE HEHE EHH EEE IIE EEE EEE IEE BEE HIE FE IEE SESE FE EFC IEIEDE IESE IIE IE SEDC FE FEE SEE BEE SE IESE DE 


* * 
* * 
* * 


PEK IEHEHE FEFE IE IE HE FE SE FE IEEE HE IE IEE IE FE IE IESE FEE FETE HEEL HEHEHE IEFEEDE DE FE IE FE IE IESE FE HE IESE FE IESE FE IESE FE IEEE IEE HEHEHE IEE 
PEK HEHE HEE FEE HEHEHE HEHE IEE ESE IE HE ESE HEINE HEE HE EEE HE IE FEE FE EIEN FEED FE SEE SESE FE FE FE HEHE HEHE EEE IESE HEE HEHE HEHE 


x % 
* GROUP 7. VARIABLES STORED, SOLVED AND NAMED * 
% * 
SR HE FE IE BE BE HE EAE IEEE FG HE HE BE OE MES SEH MERE IESE TEN IER OE SIO IIE I BE 
x w 


SOLVE(P1,H1,U1,V1,W1) 
SOLUTN(P1,Y>Y>,Y>NoNN) 
STORE( RHO1,TMP1 ) 

* 

* 

% 


* 
* 
* 


PEPE IED SESE HEHEHE SESE SESE FE FE HE IESE SE IE ESE SE HEE ESE HE ESE SE FEE HE IEE IEICE IESE IESE HE FETE HE HE FETE EEE ETE FE FE FE FEE IEE HE ETE 
HHH IA IH HEH HEHE HEHE FEE EEE IE EE HE HEHEHE IE HEE HEHE HEISE HIE HE IEFE ESE HEIEFE HEH HEE IESE HE EEE FE IE FIE IE HE HE FEE HEHE HEHE HE 


* 


* GROUP 8. TERMS (IN DIFFERENTIAL EQUATION) AND DEVICES 


* 


* 
* 
* 


FE FEE IE IE EEF HE EE FEE HEHE FEE IEE HE IEE IE EE FEE EEE IEA IEE HE EE HE HEISE HE IEE IEICE HE HEE HE ENE HEE SE HE HEHE HK 


DEFAULT TAKEN 


KK KK kK 


KK KK K 


PEI FEI IE IEE HEME HE EE FE HEE IE EE IEE HE IE IEE ESE IEE HE HE IE IEE IE EE IE HEE EE EEE HEE HE IE HE HEA IE HE IEE HEE HE EEE HEE HEHE HEHE 


% 
% 
% 


HEE IEE ETE FE HE EI HE IE HE EE IEE EIT FEE FEE EE HEI ESE HEHEHE ESE IE HE FEE TE HEE ESE SE FE HE HEE IEE EME HE HEE EEE HE HE HE HEHE 


~ 
* GROUP 9. PROPERTIES OF THE MEDIUM (OR MEDIA) 
* 
* 
TMPLA=TINY 
TMP AMB=295. 
PRSAMB=0. 
CPAMB=1004. 
TMP 1B=1.0/CPAMB 
TMP 1=GRNOD2 


TMPJET=2000. 
CPJET=1004. 
PRSJET=0. 
GASTCON=287. 
RHO1B=1./GASCON 
MASJET=100. 
BReoou-101S25. 
RHOJET=( PRESSO+PRSJET J/( GASCONXTMPJET } 
VELJET=MASJET/RHOJET 
RHO1=GRNOS 
ENUL=1.0E-05 
DRHIDP=GRNDS5 
TURMCD(KEMODL ) 


oye, 


* 


STORE (VIST ) 

KEINJ=0.5*(0.005*VELJET )x*2 

EPINJ=0. O9*KEZINJ**1.5/(0.05%0.5) 

KEINA=0.5*(0.005%1.0 Jx*2 

EPINA=0.09*KEINA¥*1.5/(0.1) 
x x 
x * 
FOIE EI IE IE IEE IE IEE IIE HE EIDE IE HE HE IEIEFEH FEHEIE HE EIEIENE IESE IE EEE IE HED FP TEE ETE IEE HEHE HE HEFL IE HEE 
HH HII IH IRIE HIE HEE IEEE IEE IEEE EHD HIE EIEIO IIE HII IEE 


* * 
* GROUP 10. INTERPHASE-TRANSFER PROCESSES AND PROPERTIES * 
* x 
HIE IE IIE HEE HEE HEHE IEE IEE HEIR IEEE IEE IE HIE EDIE IEE NEE IEE 
* DEFAULT TAKEN * 
x * 


FHI POH IEHE HE IESE HEE EE HE HEHE HE IEE HE IESE HE ESE HE EEE HE HE HIE HE HEHEHE HEIL DED HEHE IEDE HE HEE DEE HEME IEHE FEE HEHEHE HE HE 


* * 
* GROUP 11. INITIALIZATION OF VARIABLE OR POROSITY FIELDS * 
x “ 
SEK HIE IE A EI IEEE IE KE FE SESE IEEE IEE SE FE IEE IE EE IE IEE FE EI EE ge Se oh IE EEE EE HK 
* * 


FIINIT(H1 })=CPAMB*XTMPAMB 
FIINIT{ TMP1 }=TMPAMB 
FIINIT(DENI J=1.0 
FIINIT(P1)=0.0 
FIINIT(U1)=0.0 
FIINIT(V1J)=0.0 
FIINIT(W1)=0.0 
FIINIT{KE J=KEINA 
FIINIT(EP }=EPINA 

x 

+ 

x 

* BLOCKING OUT THE MASS OF THE JET ENGINE 

* 

CONPOR(0.0,CELL,1,1,5,6,6)7) 

x % 
WALL,FLOOR, AND CEILING BOUNDARY CONDITIONS. RESTRICTING FLOW* 
ACROSS A BOUNDARY BY USING CONPOR VALUES FROM 0.0 (NO FLOW) TOx 
1.0 (MAX FLON). BOUNDARY FRICTION TO FLOW IS ALSO IMPLEMENTED 
USING A NEGATIVE BOUNDARY COORDINATE. THIS CONVENTION CAUSES * 
FLON PARALLEL TO THAT BOUNDARY FACE TO HAVE FRICTION. BOUNDARY* 
CONDITIONS WILL BE DESCRIBED BY COMMENT AND PLANE. * 


KK KK * 


KK KK KK K KK 


TEST FACILITY FOUNDATION AND GROUND OUTSIDE STRUCTURE 


x KK * 


x 
CONPOR(0.0,SOUTH,1,NX51,1,1,NZ) 
PATCH( FLOR ,SWALL,1,NX,1,1,1,NZ,1,1)3 
COVAL(FLOR,U1,GRND2Z,0.03 
COVAL(FLOR,W1,GRND2,0.0)} 

COVAL( FLOR,KE »GRND2 »>GRND2 } 


COVAL( FLOR,EP »GRND2Z »GRND2Z } 
x * 


* JET-ENGINE TEST BAY 

: BAY SIDE WALL (Y-Z) 
GONPORI0.0,HEST» 95991125259) 

: BAY CEILING (xX-2Z) 
Peosione.coUrH ava, 120s.2, 5) 

: BAY WALL REAR (X-Y¥) 2 SECTIONS 

* 


CONPOR(0.0,LOW,6,8,1,12,-105,10) 
CONPOR(0.0,L0W,1,5,11,12,-10,10) 

* 

* EXHAUST DUCT INLET TUBE 

* 
CONPOR(0.0,NORTH,1,2),3,3,105,14¢) 
PATCH( TUBL,NWALL ,1,2,3,35,105,1451)1) 
COVAL(TUB1,U1,GRND2,0.0) 
COVAL( TUB1,W1,GRND2,0.0) 
COVAL( TUB1,KE ,GRND2,GRND2Z) 
COVAL(TUB1,EP,GRND2,GRND2 ) 
CCONPOR(0.0,S0UTH,1,2,8,8,10,14) 
PATCH( TUS2,SWALL»,1,25,8,8,10,14,151) 
COVAL(TUS2,JU1,GRND2,0.0) 
COVAL( TUB2,W1,GRND2,0.0) 
COVAL(TUB2 ,»KE »GRND2,GRND2 ) 
COVAL(TUB2,EP,GRND2,GRND2) 
CONPOR(0.0,WEST,3,354% 57,105,149) 
PATCHI{ TUBS »NNALL »3 5354 97,10,1451,1) 
COVAL(TUB3,V1,GRND2,0.0) 
COVAL(TUSS  ,WN1,GRND2,9.0) 
COVAL( TUB3 »>KE »GRND2 »GRND2 ) 
COVAL( TUB3,EP,GRND2,GRND2 ) 
CONPOR(G.0,LOW,1,55,1,5,10,10) 
CONPOR(0.0,LOW,1,5,8,10,10,10) 
CCHPOR(0.0,LON,2,55957,10,10) 

* 

* SECONDARY-AIR INLET WALL (X-Y) 

od 
CONPOR(0.0,LCW,6,7,1,10,12,12) 
PATCH(SAIW, LWALL,6,7,1,10,12,12,1)1) 
COVAL( SAIW,U1,GRND2 0.0) 
COVALISAIW,V1,GRND2,0.0) 
COVALISAIW,KE,GRND2,GRND2 } 
COVAL( SAIN, EP,GRND2,GRND2 ) 

* 

- SECONDARY-AIR INLET AND PARTIAL EXHAUST ODUCT CEILING (X-Z) 

x 
CONPOR(0.0,SOUTH,1,7,11,11,10,12) 
PATCH(SAIWPE ,SHALL,1,7,11,11,10,12,1,1) 
COVAL(SAIWPE ,U1,GRND2,0.0) 
COVAL(SAIWPE »W1,GRND2,0.0) 
COVAL(SAIWPE ,KE ,GRND2,GRND2Z ) 
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* 


* 


COVAL{(SAIWPE ,EP,GRND2,GRND2 ) 

% 

* EXHAUST DUCT CEILING (X-Z) 

* 
CONPOR(0.0,SOUTH,1,5,11,11,13,20) 
PATCH( EDC,SWALL»,1,5,11,11,13,20,1,1) 
COVAL( EDC ,U1L,GRND2,0.0) 

COVAL{( EDC ,W1,GRND2,0.0) 
COVAL( EDC >KE ,GRND2,GRND2 ) 
COVAL( EDC ,EP,GRND2,GRND2 ) 

* 

* EXHAUST-DUCT WALL (Y-Z) 

% 

CONPOR( 0.0 ,HEST,6,6,1,10,13,20) 
PATCH( EDNXZ,WNALL »6.651,10,13,20,1,1) 
COVAL( EDWXZ ,»V1,GRND2,0.0) 

COVAL( EDWXZ,WH1,GRND2,0.0) 

COVAL( EDWXZ KE »GRND2,GRND2 ) 

COVAL( EDWXZ,EP,GRND2,GRND2 ) 

x* 

ve EXHAUST-STACK WALLS (X-Y) THO SECTIONS 

x* 
CONPOR(0.0,HIGH,1,5,11,13,20,20) 
PATCH( ESNXY1,HNALL>1,5,11,13,20,205,1)1) 
COVAL( ESNXY1,U1,GRND2,0.0) 

COVAL( ESNXY1,V1,GRND2,0.0) 

COVAL( ESWXY1,KE ,GRND2,GRND2 ) 

COVAL( ESNXY1>,EP,GRND2,GRND2Z ) 
CGNPOR(0.0,LOWN,1,5,1,13,27,27) 
PATCH(ESNXY2,LWALL,1,5,1,13,27,27,1),1) 
COVAL( ESNXY2 ,U1,GRND2,0.0) 

COVAL( ESNXY2,V1,GRND2,0.0) 

COVAL( ESNXY2,KE »,GRND2,GRND2 ) 
COVAL(ESWNXY2,EP,GRND2,GRND2 ) 

* 

* EXHAUST-STACK WALL (Y-Z) 

pa 
CONPOR(0.0,WEST,6,6,1,13,21,26) 
PATCH( ESHXZ ,WWALL»,6,6,1,13,21,265,1)1) 
COVAL(ESHMXZ ,V1L,GRND2,0.0) 

COVAL( ESHXZ,W1,GRND2,0.0) 
COVAL(ESNXZ,KE »GRND2,GRND2 ) 
COVAL( ESHXZ,EP,GRND2»,GRND2 ) 

* 

* TURNING VANES 

7 
CONPOR(0.0,SOUTH,1,55-7,7,21,22) 
CONPOR(0.9,SOUTH ,1,5,-3,3,23 5,24) 
CONPOR{ 0.0,LON,1 Fost s lOb=25 »23 ) 
COHNPOR(0.0,LOW,1,5,3,6,-25,25) 


*K 


* 


*K 


* 


* 
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# GROUP 12. CONVECTION AND DIFFUSION ADJUSTMENTS 


* 


38) 


x 
* 
x 
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* DEFAULT TAKEN 
“ 


* 
% 
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* GROUP 13. BOUNDARY CONDITIONS AND SPECIAL SOURCES 


* 


* 
* 
% 


HHH IIH HHI IIE HE HEH HHH HK HK IIH HHH IH IIE HE HE HE EE IEE HIE IEE IEEE FEE IE IE IE HE EEE FEES IE 


x 
* JET-ENGINE MASS SOURCE 
x 
PATCH( JETIN,HIGH,1,1,5,6,5,5)1)1) 
COVAL(JETIN,P1,FIXFLU,-MASJET*1. ) 
COVAL( JETIN>,H1,ONLYMS »SAME ) 
PATCH( JETOUT ,LOW,1,1,5,6,8,8 5151) 
COVAL( JETOUT ,P1,FIXFLU,MASJET*1. ) 
COVAL( JETOUT ,W1,ONLYMS ,VELJET ) 
COVAL{( JETOUT »H1,ONLYMS ,;CPUETXTMPJET ) 
COVAL( JETOUT ,KE ,ONLYHS »KEINAJ } 
COVAL( JETOUT ,EP ,ONLYMS ,EPINJ ) 
* 
ae ATMOSPHERIC BOUNDARY CONDITIONS 
x 
* LAST CELL SURFACE IN X DIRECTICN 
*% 
PATCH( XSKY , EAST »NX NX 1 NY ,15,NZ,1>513 
COVAL( XSKY »P1,0.1,0.0} 
COVAL( XSKY >H1,ONLYMS ,CPAMB*TMPAMB ) 
COVAL( XSKY »KE ,ONLYMS »KEINA ) 
COVAL( XSKY »EP,ONLYMS,EPINA ) 
* 
* LAST CELL SURFACE IN THE Y DIRECTION 
x 
PATCH( YSKY »NORTH,1>NX NY ONY, 1 NZ,1)51) 
COVAL( YSKY »P1,0.1,0.0) 
COVAL( YSKY,H1,ONLYMS » CPAMB*TMPAMB ) 
COVAL( YSKY,KE ,ONLYMS ,KEINA ) 
COVAL{ YSKY,EP ,ONLYMS ,EPINA ) 
* 
* FIRST CELL SURFACE IN THE Z DIRECTION 
* 
PATCH({ ZSKY1,LOW,1 NX, ,L SNY,1,1,1,1) 
COVAL( ZSKY1,P1,0.1,0.0) 
COVAL( ZSKY1,H1,ONLYMS ,CPAMS*XTMPAMB ) 
COVALI ZSKY1,KE ,ONLYMS ,KEINA } 
COVAL( ZSKY1,EP ,ONLYMS ,EPINA ) 
x 
* LAST CELL SURFACE IN THE Z DIRECTION 
x* 
PATCH ZSKY2,HIGH,1,NX 51 ,NY »NZ»NZ>,1513 
COVAL( ZSKY2,P1,0.1,0.0) 
COVAL{ ZSKY2,H1,ONLYMS ,CPAM3*TMPAMB ) 
COVAL( ZSKY2 > KE »ONLYMS »>KEINA ) 


ayy 


% 
% 
* 
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COVAL(ZSKY2 ,EP ,ONLYMS ,EPINA } 
SE SEM ot ae EE EEE A HEHE EH HEHE EH IEEE A HEE IEE BE HE IIE DE FEI 3 IEIE FE He IEEE BEE IE IEICE 


FEE IE HE DE FE HEHE TE FE HE FEE FEE IEE QE HEHE FE FE FEE TE FEE FE HE EEE HE FEE EEE FEE FEE FE FE FETE EEE SHE FEES FE HEE HE FEE EEE HE HEHE FE 


x * 
* GROUP 14. DOWNSTREAM PRESSURE FOR PARAB=T * 
* 2 
SMR EH EI HIE KHONG SEI RIES HI IOE IIE IEE IES III OI GEE IR a oe 
* ” 
* DEFAULT TAKEN * 
~ “ 
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HEE IE ETE HE HEFIE HE HEHE EEE HEH EEE HEHEHE HEHEHE IEE IE FETE HE FE FE FE HE HE FE IE FE FE FE FE HE IEE HE EE FEE IESE FEE IESE FE IESE EH HEIHEIOIE NE 


% * 
¥ GROUP 15. TERMINATION OF SWEEPS * 
% % 
YEHE FEE FEN FE HE EH IEE EEE HEE HEE THE EE EEE HEE ETE FEE E HEHEHE IEEE FOIE FE FEE FEE FEE FE FEE FE FE FE FE EEE FEE EEE HE FEE IESE 
* * 

RESTRT(ALL ) 

FSWEEP=1 

LSWEEP=2 
% * 
% * 
% * 


MH IE IE HE IE HE HEHEHE HEHE FEE HE HEHE HEE FE FEE IEEE HEE EEE HE HE IEE IE HE HEE EE HEHE IEF FE HEHE E FE HE HEHE HEHEHE HEHEHE HEHE HEHE HHH 
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% % 
% GROUP 16. TERMINATION OF ITERATIONS * 
* * 
KH IEIE FEE HE IEE FE HEE IE EEE HE HEE SE HE IEE FEE FE FE HE HEE IEE FEE IE FE FE FE FE HE IE TE HE FEE EEE FEE EEE FE FESE TEE FEFEDSE HE FE HE HEHEHE 
x % 
% DEFAULT TAKEN * 
% * 


FEE FE HEHEHE HE FETE FEE HE FETE HE IEE FETE HE FIE HE HE FETE FETE HE HE FETE FEE HE FE HE HE IEE HE FE FE ETE FE FEE EI HIE HE FE FEE HE HE FE SE HEHE FEE HE HEHEHE 
MEH HK HEHEHE HEHE FE IE HEH FEE HE ETE HE HEHE FETE HE HE IEE IEE HEHE HE FE FETE HE FE FETE HE ETE EE FE FE FE HE FE FE FE HE ESE HE HEEL IE IE IE FE FE HE FEE 


x * 
* GROUP 17. UNDERRELAXATION DEVICES % 
* * 
IOI HE IIE EH EE IE FEF AE IE SE EEE HEFL ESE EE FE IESE IE ESE HIE IE FOIE HEHE FETE EFC DEE IEE IK HE FE IIE HE 
~ * 
DELT=50. /(NZ*0. C¥VELJET } 
FACT=0.25 
SHIFTED TO FACT=0.5 AT 1550 SWEEPS 
FACT=0.5 


RELAX(P1,LINRLX,0.3 ) 

RELAX(U1,FALSOT,FACT*DELT } 
RELAX(V1,FALSOT,FACTXDELT } 

RELAX(W1,FALSDT ,FACT*DELT } 

RELAXUKE ,FALSDT »FACTXDELT } 

RELAX(EP,FALSOT ,FACTXDELT } 

REMOVED THE 2.5 FROM H1 RELAX AT 5000 SWEEPS 

RELAX(H1,FALSOT ,FACT*DELT¥1.03} 


* * 
%* * 
* * 
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o * 
* GROUP 18. LIMITS ON VARIABLES OR INCREMENTS TO THEM * 
* a 
HEHE REIEIE HE HEHE HEI IE HEHEHE HE HEHEHE HEHE HEHEHE IEEE HIE IIIA EI HHI IIE 
“ % 


VARMIN(P1 )=-2*PRESSO 
VARMIN(H1 )=CPAMB*TMPAMB 
VARMIN( U1 )=-VELJET 
VARMINCV1 J=-VELJET 
VARMINOW1 J=-VELJET 
VARMIN( TMP 1 )=TMPAMB 
VARMIN(€ RHO1)=0.1 
VARMAX(P1 )=2*PRESSO 
VARMAX(CH1 J=CPJE TXTMPJET 
VARMAX( U1 J)=VELJET 
VARMAXCV1 J=VELJET 
VARMAX(W1 J=VELJET 
VARMAX( TMP1 )=TMPJET 
VARMAX(RHO1 J=2.0 


* * 
* * 
* * 
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* * 
* GROUP 19. DATA COMMUNICATED BY SATELLITE TO GROUND * 
* 4 
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x “ 
% DEFAULT TAKEN * 
* * 
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* GROUP 20. PRELIMINARY PRINTOUT * 
~ « 
HH HHH HHH HAH IH SHAIKH HAHAHAHAHA IESE SISSIES HAH IH 
* * 
~ 
* DEFAULT TAKEN * 
* * 
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* GROUP 21. PRINTOUT OF VARIABLES * 
* # 
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* * 
* DEFAULT TAKEN * 
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* GROUP 22. SPOT-VALUE PRINTOUT * 


cn 
(D 


IXMON=1 
IYMON=5 
IZMON=22 
TSTSWP=10 

%* 

x 

x 


* 
* 
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* GROUP 23. FIELD PRINTOUT AND PLOT CONTROL 


* 


* 
* 
* 
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* 
PATCH( PLOTS,PROFIL,1,1,5,55,1,NZ ,1,1) 
PLOT(PLOT5,P1,-3000.0,2500.0) 
PATCH( PLOT3,PROFIL,1,1,5,5,1,NZ51,1) 
PLOT( PLOT3 ,»U1,-20.0,50.0) 
PATCH( PLOT4,PROFIL,15,1,5,5,1,NZ5151) 
PLOT(PLOT¢,V1,-40.0,60.0) 
PATCH( PLOT2 »PROFIL,1,1;5,5,1,NZ,1>)1) 
PLOT(PLOT2Z>,W1,-15.0,500.0) 
PATCH(PLOT1,PROFIL,1,1,5,5,1,NZ,1)51) 
PLOT(PLOT1,H1,275000.0,2100000.0) 
PATCH( PLOT6,PROFIL,1,1,1,10,20,20,1,1) 
PLOT(PLOT6,H1,300000.0,1000000.0) 
PATCH( PLOT7,PROFIL,1,1,1,10,20,20,1,1) 
PLOT(PLOT7>,A1,20.0,175.0) 
PATCH( PLOT8,PROFIL,1,5,5,5,20,205,1,1) 
PLOT( PLOT8,H1,300000.0,1000000.0) 
PATCH( PLOT9,PROFIL,1,5,5,5,20,20,1,1) 
PLOT(PLOT9,W1,20.0,175.0) 
PATCH( PLOT10,PROFIL»)15,1,1,10,16,16,151) 
PLOT(PLCT10,H1,300000.0,1000000.0) 
PATCHtPLOT1I1,PROFIL,1,1,1,10,16,16,1,1) 
PLOT(PLOT11,N1,20.0,175.0) 
PATCH{ PLOT12 »PROFIL,1,5,5,5,16,165,1)1) 
PLOT(PLOT12,H1,300000.0,1000000.0) 
PATCH(PLOT13,PROFIL,1,5,5,5,16,1651)1) 
PLOT{ PLOT13,W1,20.0,175.0) 
PATCH{ PLOT14,PROFIL,1,1,13,13,21,265,151) 
PLOT(PLOT14¢,H1,500000.0,600000.0) 
PATCH{ PLOT15,PROFIL>»1,1,135,13,21,2651,1) 
PLOT(PLOT15,V1,10.0,45.0) 
PATCH( PLOT16,PROFIL»,1,5,13,13 5235235151) 
PLOT(PLOT16,H1,500000.0,600000.0) 
PATCH(PLOT17,PROFIL,1,5,13,13,23,235151) 
PLOT(PLOT17,V1,10.0,45.0) 
NUMCLS=1O*N*PRIN 
NPLT=10 
ITABL=3 
IPROF =3 


LUPR3=6 
“ 
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* GROUP 24. DUMPS FOR RESTARTS 
x 


* 
* 
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Figure C.2. Staggered-Grid Cell Configuration 
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Cartesian Cell Configuration 


30) 


HIGH 


RUN TITLE AND OTHER PRELIMINARIES 

TE SoLENCE TIMeE-STEPASPECIFICATIONS 
X-DIRECTION GRID SPECIFICATION 
TeDoRECtlON@Grap SPECIFICATIONS 
ZeolRectlONeGR ID SPECIFIGATIONS 

DOU er lIreEDeCOURDINATES GREGKID DISTORTION 
VARIABLES STORED, SOLVED AND NAMED 

TERMS CIN DIFFERENTIAL EQUATIONS) AND DEVICES 
PRORPERI IES SOR THE MEDIUM (CGR MEDEA) 
PNTEREHeASE=>RANSEER FROGESSES AND PROPERTIES 
INITIARI ZATIONS Or VARTABEEROR POROSITY FIELDS 
CONVECTION AND DIFFUSION ADJUSTMENTS 
BOUNDARY CONDITIONS AND SPECIAL SOURCES 
DOWNSTREAM PRESSURE FOR PARAB=T 

TERMINATION OF SWEEPS 

TERMINATION OF ITERATIONS 

UNDERRELAXATION DEVICES 

CIMITSmGN VARIABLES OR INGREMENTS 10 THEM 
DATA COMMUNICATED BY SATELLITE TO GROUND 
PRELIMINARY PRINTOUT 

PRINTOUT OF VARIABLES 

SECTSVALUC PRN hour 

FIELD PRINTOUT AND PLOT CONTROL 

DUMPS FOR RESTARTS 


Figure C.4. PHOENICS Data Group Catagories 
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Bie: COPYQ@1] EXE Al 


&TRACE 
Solar l 
&TYPE WHAT ARE THE FILENAME AND FILETYPE OF THE Q1 FILE YOU WANT 
capcoee TO STORE O88 MVS? USE ENTER TO EXIT. 
&READ VAR &Q1 &Q2 
SIF 7&Q1 EQ 7 &EXIT 
STATE &Q1 &Q2 A 
SIF &RC EQ 0 &GOTO -NEXT 
Abpe sr lUEesdlesQez Ae DOES HON EXIST. USE EXTER TO EXIT. 
&READ VAR & 
SIF 7&8 EQ 7 &EXIT 
&GOTO -START 
=HEXT 
STATE COPYQ1 JCL A 
Sif &hC seGe Oe SeoloOy-COPyT 
TYRE Teele PlEE CORT GQ! JCL A DOES NOT EXIST ; 
&STYPE SEE LAB MANAGER. 
SEXI Taae 
meOr Yy 
SERROR &EXIT -1 
Coir Ite COPYQI JCE A TEMP1 JCL A CREPL 
& = &CONCAT OF /ZYOURQIFILENAMNE/Z &Q]1 7 
&STACK TOP 
Sone k CLOCATE7YOURQIFILENAME7 
&STACK CHANGE & 1 1 
Sornck CLOCATE/SYSUTI~ 
&SSTACK GET &Q1 &Q2 A 
&STACK TOP 
&STACK FILE 
FEDITeTEnPL JCL A 
Exec SUBMIT TEMP1 JCL A 
CONHWAIT 
DESBUF 
S&S EXIT 


meee: COPYQ1 JEL Al 


ZACANUSH JOB €49555,9999), 'COPYTMVS',CLASS=A 

Zee EXEC PGM=ITEBGENER 

77SVSPRINT DD SYSOUT=A 

ZeovoiN DD DUMMY 

meer ounce DD DSH=MSS .S4555.PHOENITCS.SRECYOURQIFILENAME),DISP=SHR 
//SYSUT1 DD X 

oes 

17 


Figure C.6 COPYQ]1 EXEC And JCL Files 


FILEs RUNS AT EXEC Al 


& TRACE 
-START 
&TYPE WHAT IS THE FILENAME OF THE Q1 FILE YOU WANT TO RUN CSATELITE)? 
&TYPE USE ENTER TOME cei 
&READ VAR &Q] 
&IF 7&Q] EQ A SExIT 
* STATE &Q] DATA A 
* &IF &RC EQ 0 &GOTO -NEXT 
* &TYPE FILE &@1 DATA A DOES HOT EXITS. USE EHIER [082 
4 &READ VAR & 
* &IF 7& EQ / &EXIT 
* &GOTO -START 
-HEXT 
STATE RUNSAT BCH A 
&IF &RC EQ 0 &GOTO -COPY 
&TYPE TEMPLATE FILE RUNSAT JCL A DOES TUT EX est y; 
&TYPE SEE LAB MANAGER. 
SEX lie 2 
-CUrY 
&ERROR &EXIT -1l 
COPYFILE RUNSAT JCI A TEMPZ JEL A CREPE 
& = &CONCAT OF 7 YOURQIFILENAME 7 &Q] / 
&STACK TOP 
&STACK CHANGE & * Il 
&STACK FILE 
XEDIT TEMRaee Cie 
EXEG SUBMI I BREMEze Cl FA 


CONWAIT 
DES BUE 
SEXIT 


Figure C.7”~ RUNSAT) EXE@ Frigre 


34 


FILE: RUNSAT AGE Al 


//7RXHUSH JOB (4555,9999), 'XCHAM', CLASS=A,REGION=1024K 
7/GO2 EXEC PGM=SATLITP 

//STEPLIB DD DSN=MSS.PHOENICS.LOAD, DISP=SHR 
faut sti DD DSN=MSSaRHOENTES UIST. DIR, DISP=C(0 
Zee) 2 DD DSN=NSss: PHORBMICS.CIST2. DIR, DISP=C(O0 
meeA STS DD DSH-MSS .PHOENICS LIST3.DIR,DISP=C(0 
//LISTG DD DSN=MSS.PHOENICS.LIST4.DIR, DISP=COLD 
//Q1 DD DSN=MSS.S4555.PHOENICS .SRCCYOURQIFILENAME),DISP=COLD) 
//PHHELP DD DSH=MSS.PHOENICS.PHHELP.DATA, DISP=COLD) 


LD) 
LD) 
LD) 
) 


MAOCRL DD DSN=MSS.PHOENICS.SCRL,DISP=CNEW, PASS), 
as DCB=(RECFM=F,LRECL=80, BLKSIZE=80), 

// SPACE=S(TRK,C12,3)),UNIT=SYSDA 

ee COry DD DSN=MSS.PHOENICS.COPY, DISP=(NEW,PASS), 
// DCB=(CRECFM=F,LRECL =80,BLKSIZE=6520), 

a7 SPACE=(TRK,(€12,3)),UNIT=SYSDA 


Zan O9FO0T DD sy souT=x 
//FTO6FO001 DD SYSOUT=x 
//FT14F001 DD SYSOUT=*x 


76512 
LAGOS 
4/6514 
775515 
4/6516 
y7G6518 
4/6526 
4/6520 
76521 
7 OSCe 
WeGoe 3 
4/6524 
4/6527 
//6528 
4/6529 
4/6530 
4/6531 
eG) 55 
776156 
ZAG) S37 


DSN=MSS. 
DSN=MSS. 
DSH=MSS. 
DSH=MSS. 
DSH=MSS. 
DSN=MSS. 
DSH=MSS. 
DSH=MSS. 
DSN=MSS. 
DSH=MSS. 
DSN=MSS. 
DSN=MSS. 
DSN=MSS 
DSN=MSS. 
DSH=EMSS. 
DSNEMSS. 
DSH=MSS. 
DSN=MSS. 
DSH=MSS. 
DSH=MSS. 


PROENTCS. 
sO 135, DIS) 


PHOENICS 


0 
0 
Pa SCS; 0 
PHOENMICS. 0 
iphel aie Gs, 0 
RmOeiIrc >. 0 
PaO Chiles. 0 
-6520,DISF-c0 
~G521,08SE-<0 
:6522, DIsk=e0 
Q 
0 
0 
0 
0 
0 
0 
0 
0 


BnOENIECS 
PHOENICS 
PROENICS 


PHOENTCS: 
PHOENTCS: 
sRHOEMNLES G52/, Dist 
ipl 8; alc oes, & 
PHOEMICS: 
PHOENITCS< 
PHOENICS: 
EmOEMIC Se 
PROENTCS:. 
ec7o/ » DiSoE=(e 


PROENILCS 


Ciic, US io. 


G514, DISP=( 
Gol 5, Dlsra« 
GS16,DiSr =< 
Gale, Urs =« 
G526,DISP=¢ 


G523,DISP=¢ 
G524¢,DISP=( 


G528,DISP=( 
G529,DISP=¢ 
G530,DISP=C 
Goo, DIS = 
G735,DI3P=(¢ 
Grso7DISk=¢ 


-PROENTCGs ~DEuG 
-FHOENT GS phe sD SP=SHR 


DD DSN=MSS.$4555 
DISP=SHR 
DD DSN=MSS .34555 


S/DF10 
a 
yer Le 
ae 


Figure C.8. @sRUNSAY JCL File 
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FILE: RUNEAR 


& TRACE 


EXEC Al 


EXEC SUBMIT RUNEAR JCL A 


FILE: 


77 EXHUSE 


RUNEAR 


wel Al 


JOB (4555,9999), 'CHAM', CLASS=G, REGION=2048K 


//¥MAIN ORGSNPGVM1] .4555P,LINES=(C12) 
//GO2 EXEC PGM=EARTHP 
//STEPLIB DD DSIN=11SS.PHOEWICS LOAD ROT eon 


A/A DS TL DD DUMMY 

//DFOYI DD DSH=MSS.54555.PHOENICS.DF09, DISP=SHR 
77 VEOL DD DSN=SSDFE01 2 DISP=CHEWN FP ASss>: 

TZ DCB=(RECFM=VBS,LRECL=19000,BLKSIZE=19006), 
V7, SPACE=CIRK 62523 )) UIT -s 50-4 

AADEDZ DD DSN=&&DFO2, DISP=CHEHM, PASS), 

es és DCB=(RECFM=VBS,LRECL=19000,BLKSIZE=19006), 
TA SPACE=(TRK. C25,-5)) > UHM =s ona 

//DFO3 DD DSH=&&DF03,DISP=CNEM, PASS), 

We DCB=(RECFM=VBS,LRECL=19000,BLKSIZE=19006), 
Vite SPACE=CTRK,€25, 3)),UNTheSYSDA 

//DFEOG DD DSN=&&DFO0G, DISP=(NEW, PASS), 

Ve DCB=(RECFM=VBS,LRECL=19000, BLKSIZE=19006), 
// SPACE=(TRK,(25,3)),UNIT=SYSDA 

77 DFOS5 DD DSN=&&DF05,DISP=C(HEMW,PASS), 

Yat DCB=(RECFM=VBS,LRECL=19000, BLKSIZE=19006), 
// SPACE=(TRK,(€25,3)),UNIT=SYSDA 

//DFEO7 DD DSN=&&DF07,DISP=(NEW,PASS), 

TT DCB=(RECFM=VBS,LRECL=19000, BLKSIZE=19006), 
Bas SPACE=(TRK,( 25, 3)), Uli -sSyYs0e 

//DFO8 DD DSH=&aDFE03, DISF-CHEN.  PAsos 

VA A DCB=(RECFM=VBS,LRECL=19000, BLKSIZE=19006), 
SAS SPACE=(TRK,(25,3)),UNIT=SYSDA 

77 DENS DD DSN=&&DF15,DISP=CHEW, PASS), 

VA A DCB=(RECFM=VBS,LRECL=19000, BLKSIZE=19006), 
Ve SPACE=(TRK; (25,30), 0UN1I-=SvsvA 

//DF16 DD DSH=MSS.S84555.PHOENICS.DF16,DISP=SHR 
//DF10 DD DSN=MSS.84555.PHOENICS.DF10,DISP=COLD) 
//DFi2 DD DSN=MSS.846555.PHOENICS .DFIZ, Disr-(OLlbD 
//FT1GFOO01 ODD SYSOUT=*x 

//FTO6FO001 DD SYSOUT=x 

tf, 


Figure C.9 RUNEAR EXEC And JCL Files 
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Test Facility Side View (Y-Z Plane) 
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Test Facility Top View (X-Z Plane) 
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Figure C.12. Vertical Y-Z Plane Zones 
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Figure C.13. Horizontal X-Z Plane Zones 
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* GROUP 3. X-DIRECTION GRID SPECIFICATION 

* NX = THE NUMBER OF CELLS IN THE X DIRECTION CINTEGER) 
* XULAST = DISTANCE IN THE X DOMAIN (CMETERS) 

X 


NX=XNCBL1+XNCMD1+XNCMD2+XNCMD3+XNCBL2+XNCBL3+XNCEX) 
XULAST=XLBL1+XLMIDI+XLMID2+XLMID3+XLBL2+XLBL3+XL EXT) 
® 


XFRACC1 )=-XNCBL1; XFRACC2)=XLBLI/“CXNCBLIXXULAST) 
XFRACC3)=XNCMDI; XFRACCG)=XLMIDI/ CXNCMDIXXULAST ) 
XFRAC(C5) =XNCMD2; XFRAC(C6)=XLMID2/ CXNCMD2¥XULAST) 
XFRACC7)=XNCMD3; XFRACC38)=XLMID3/ CXNCMD3XXULAST ) 
XFRACC9)=XNCBL2; XFRAC(C10)=XLBL2/7 CXNCBL2¥*XULAST ) 
XFRACC11)=XNCBL3; XFRACC12)=XLBL3/ CXNCBLSXXULAST ) 
XFRACC13)=XNCEX1; XFRACC1G)=XLEXT1/Z CXNCEXIXXULAST) 
x 


* GROUP 4. Y~DIRECTION GRID SPECIFICATION 
x NY =SIRe NUMBER SOR CELE S IN THE YT DIRECTION CINTEGER) 
*X YVLAST = DISTANCE IN THE Y DOMAIN CMETERS) 


x 
nV=YHCMID+YNCBLI+YNCBL2+YNCEXI+YNCEX2+YNCBLS 
feel -YLMID+YLBLI+YLBL2+YL EXTIFYLEXT2+YL BLS 

® 


YFRAC(C1)=-YNCBLI; YFRACC2)=YLBLI/SCYNCBLIXYVLAST ) 

YFRACC3)=YNCMID; YFRACCG)=YLMID/CYNCMI DXYVLAST ) 

WmeRACCS)=YNCBL2; YFRACC6)=YLBL2/(CYNCBL2XYVLAST ) 

YFRACC?/ J)=YNCEX1; YFRACC8 )HYLEXTIZCYNCEX1¥YVLAST) 

YFRACCS )=YNCEX2; YFRACC1O)=YLEXT2/7 (CYNCEX2*®YVLAST ) 

fermacCl1)=YNCBL3; YFRACC12)=YLBL3/CYNCBL3XYVLAST) 
x 


SOGROUP SS. Z-UIKECTION CR RUeopeeteen | LON 
ex NZ = THE NUMBER OF CELE SIN Tee eZeDIRECTION CINTEGER) 
* ZWLAST = DISTANCE IN THE Z DOMAIN (METERS) 


x% 
NZ=ZNCEX1L+ZNCINI+ZNCIN2+ZNCIN3+ZNCING+ZNCIN5S+Z2NCIN6+ZNCEX2 
Maeno l=ZLEXTL+ZLINTL+ZUINT 2+Z0 INT S+ZLINGG+Zi INT S+ZLINT6+ZLEXT2 


x 
ZFRACC1)=-ZNCEXI]; ZFRACC2)=ZLEXTISCZNCEXIXZWLAST) 
ZERAC(3)=ZNCINI; ZFRACCGIEZLINTIZCZNCINIXZWLAST ) 
ZFRACC5)=ZNCIN2, ZFRACC6IVEFZLINT2“CZNCINZXZWLAST ) 
ZFRAC C7 DFZHCING; ZFRAC(8)=ZLINT3/CZNCINSXZWLAST ) 
ZFRAC(C9)=ZNCING; ZFRACCLOV=EZLINTG/CZNCINGXZWLAST ) 
ZFRACC11)=ZNCINS; ZFRACCI2)=ZLINTS/CZNCINSXZNLAST) 
ZFRACC13)=2ZNCING; ZFRACCLG)=AZLINTE6ZCZNCINGXZWLAST )D 
ZFRAC(C15)=ZNCEX2; ZFRAC(C16)=ZLEXT2/ CZNCEX2XZWLAST ) 


Figure C.14. Ql MEthod-Of-Pairs Coding Sequence 
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Let tS 2 tt 22,25, 4,5,2,2,5,5,5,5, 5,5, 5,5, 5,5, 5, 5,5, 5,5, 5 5,5, 5,555 Fatt tatotatatatetatatetetetetetetetate: 


GROUP 3. X-DIRECTION GRID SPACING 
CARTES = ik 
NX = 2 
XULAST = 9.900E+00 

METHOD OF PAIKS USED VRORSGRIO Sai inG 
XFRAC € 1) = =].Q00E+00 “GXERAC (502) = SoeGsteso2 
XFRAC € 3) = 1.Q00E+0U 5) XERAC (CG) = SS eee 
XFRAC € 5) = 1.QQQ0E+00 ;XFRAC € 6) = 5.051E-02 
XFRAC € 7) = J]. Q00E+00 GAERAC ©) g) =meS- 05 lez02 
XFRAC € 9) = 1.QO00E+00 ;XFRAC € 10) = 5.051E-02 
XFRAC € 11) = 3.Q00E+00 ;XFRAC ( 12) = 1.818E-01 
XFRAC € 13) = J OQUE+U0° | AFRAC © 15) S=ee UU 20eE-01 


HEHE HEH HE HEH HEHE HEHE HH HEHEHE HEHE HC HE HED HH HHH HEH HEH HEH DEH HHH DEH DH DH DH HD HI HE 
GROUP 4. Y-DIRECTION GRID SPACING 


NY = 16 
YVCAST = 1. 7008200 
METHOD OF PAIRS USED FOR GRID SETTING 

YFRAC € 1) = =4. 00 0E+00M WERAC Glee) = een oe ac 
YFRAC € 3) = 2.Q000E+00 ;YFRAC € 4) = 2.941E-02 
YFRAC € 5) = “G O0CE+OG 3 WERAC (lie)! = 2] oe e-02 
YFRAC © 7) = “25000E+0G = cERAC (lies) = Slee oe 
YFRAC € 9) =]. DO0E+00 4; TFRAC RGD = Se cocea. 
YFRAC € YI) = S2000E+00 | yt ERAC NC =e eo a. 


PEPE PESESE SSE SS SSS SES ESS SE SESE SESE SES SSS SESS ES SESE SS SESE SSS fF 
GROUP’ 5. Z-DIREC) ION GREDSsEACING 


PARAB = F 
NZ = 28 
ZNLAST = 6.490E+01 

METHOD OF PAIRS USED FORRGRID Seine 
ZFRAC € J) = =1.Q0Q00E+00 [ZFRAC CC” 2) = See2Ze-02 
ZFRAC € 3) =F 5.000E+00 GZERAC € 6) = 7 eeee=e2 
ZFRAC € 5) = 1.000Et00 [A2RRAC © 6) = GoaeE=iZ 
ZFRAC € 7) =3 2. 000E+00 G2ERAC Cis) = Jae 
ZFRAC € 9) == 6 .QQ0EtOU CZERAC Calo) = “ee oesi2 
ZFRAC ( TI) = 5.0Q0EF00 54 ZFRAC (32) = (Seu Zaeai2 
ZFRAC € 13) =) 6.000E+00 (|ZERAC (els) = 2eeeeeuc 
ZFRAC € 15) = 2.Q0Q00E+00 7 ZERAC Cole) = Seugzesi2Z 
ee eS SS 22S 22 22 2,2 2S 2S 22 25,22, 2,2, 2,2,0,5,5,5,5,5, 5,5, 0,5,5,5,0,2,5, 5,5, 5,5, 5,5, 5,5, 5,955,520 )3 


Pa otiieceO. | 5. 


EARTH Method-Of—Pairs Output 


ae 


1@2.0m 
3@1.8m 
5@0.5m 


4 ot et 
rr ee ee | 
—_— i ie ees ESTES 


202 .0m 


6@1.4m 


5@3.26m 
\ 


\ 

I 
; 

| 


200.75m 6@1.0m 


1@3.2m. \ 


1@2m 9@4.7m 
\ 
eae) = 
eee Pan 
| Eee 


as a a ee ae 


~< 


Figure C.16. System Domain Grid (X-Z Plane) 
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Figure C.17. System Domainecuid y-2erlane) 
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»PROEIL, 





1) 


PATCHCPLOT1O 1, , 1, 10, 16, 16, 1, 
Meese LOTI0O ,Hi , 3.000E+05, 1]1.000E+06) 
TABULATION OF ABSCISSA AND ORDINATES... 
if Hl 
Arouve-O01 3.956E+05 
fe o00=E-01 4.15GEO5 
ieecoueto0 5.46) 1E+05 
Peover w0 7.756E+05 
gaegovurOo 9,.65)0E405 
meoueErOo 9. 5S99EFOS 
feecourto0 7.667805 
Sev ouet0O 5.340/E+05 
G4.250E+00 4.097E+05 
Serovuer0O 3.613Et+05 
VARIABLE ale 
MINVAL= 3.000E+05 
MAXVAL= 1.000E+06 
CELLAV= 6.074E+05 
iid +. Parks wi ty oat sie ee en . cura oot ates re 
1 Sa se 
eo 0 + a 
ey 0 + 38 
eo 0) + a 
iO + v 
0.40 + + 
ieoo + + 
iez0 + us 
0.10 H H 
0200 Pee. th. LE cel eee oe eee: Tae he ea 
0 a me a) 4 ae. .6 ad .8 ng ~0 
THE ABSCISSA IS v6 MIN= 2.50E-01 MAX= 4.75£+00 


H vs Y-Axis at IY=1]1 to IY=10 
For IX=l1 And IZ=16 


Breure C.18. 
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1? 


PATCHC(PLOTI1 > PRORTL> (au tcnn) orl -meeel amen 1/0) imme Comune mnLEs 
PLOTCPLOT1I ,Wl , 2.000E+01, 1.750E+02) 
TABULATION OF ABSCISSA AND ORDINATES... 
Y WI 
2) S00E=-0 Len o2 SE son 
7 SOUE-01 6 Gnas! 
1 250E800 w7ale7eso. 
1 7 50c en mle eecetO2 
2,250E+00 1.696E+02 
2) 7 50 eee GET ne 
3,250E+00 1.203F+02 
3,750E+00 6.8735E+01 
G.250E+00 4.100E+01 
4.750E+00 2.919E+01 
VARIABLE W1 
MINVAL= 2.000E+01 
MAXVAL= 1.750E+02 
CELLAV= 8.653E+01 
100) hoe ee i 
OF 90s + 
0.30 + + 
Nagn = + 
60m . 
0.50 + 7 
p) ome P 
0.30 + i 
0. 200 - 
0.10 W W 
0.00 £0 a ee te ee 
0 .|l 62 63 .@ «5 .6  .7 (ee eeOMeEED 
THE ABSCISSA IS Y MIN= 2.50E-01 MAX= 4.75E+00 


W vs Y-Axis At IY=1 To IY=10 
For IX=] And IZ=16 


Froure 16..19. 


eG 


MemencPLOTIZ ,PROFIL, l, 5, 5» Dymo, Lo, l, 
PGC rPLOT12 ,Hi , 3.000E+05, 1. 000E+06 


TABULATION OF ABSCISSA AND ORDINATES... 
x 





Hl 
AeoevE-01 9.650E+05 
feeeeOr-0l 7.796Et05 
eoOer 00 5.784E405 
1.750E+00 4.848E+05 
22250700 4% .550ET05 
VARIABLE Hl 
MINVAL= 3.000E+05 
MAXVAL= 1.000E+06 
BeeeAV= 6.512E705 
P00 Ret, here hs + a ot 6 eect ie: @ ce veel rms, > Came 
te 70 t a 
teed + as 
ue + 2 
oO + at 
ee + " 
eed + + 
bao + + 
ie) + 
eo + aE 
MMP ce yt. Sy tes oot ee Cee. so Pa te ly OF 
0 I 2 3 4 , 6 7 8 -9) Lae 


Te ABSCISSA IS Xx ', MIN= 2.50E-01 MAX= 2.25E+00 
KEK KKH HK HK KKK HK KH HK KKK KKK KKK KKK KKK KKK HK KK KKK KKKKKKKKKKKKKKKKKKR 


Figure C.20. H vs X-Axis At IX=1 To IX=5 
For lY=5 And IZ=16 


Cees 


PATCHCPLOT13 ,PROFIL, lly: 5; 5; 5, 16, 16, 1, 
PLOT CEEGhIS ieee OUOPTO1 > 7 SUEF ie 


TABULATION OF ABSCISSA AND ORDINATES... 
WI 





xX 
2 7 50GE-0) Neeser. 02 
7~S0UE=0)~ lc 196E+02 
1.250E+00 6.867Et01 
1 7 S0C+007 TS. 6222500 
2.250E+00 2.5) 5250 
VARIABLE Wi 
MiNVeae=  2cOU0ETOL 
MAXVAL= 1.750Et02 
CELLAV= 8.3846E+01 
io ate Ae 
0.90 + a 
OF 0m + 
Oy ie + 
0260s + 
0.504 oF 
0.90 + + 
0. 50m + 
0.203 + 
Oe LO) oe - 
O00 Re ee ae ae I ae Me ee ee ee A G4 1 A STONE 
el 2 a 4 5 26 me no 79 
THE Sab Soe loa > X . MIN= 2.50E-01 MAX= 2225E700 


Figure C.21. W vs X-Axis At IX=1 To Ix=5 
For IY=5 And IZ=16 


HH HH HHH HH HH HHH HK HH HHH HH HH HH HHH HHH HHH HHH HH HH HH HK KH HH HK HK KKK 
PATCHCPLOT6 »PROFIL, l, l, ie; <0, 20; l, 1) 
PEO CPLOTS6 »Hl 5 SeUVUET 0S ee DOU EF0G) 


ge UN OF ABSCISSA AND ORDINATES... 


Hl 
ZeSvOE-01 5.117E+05 
meoUUE-O1 5.5356¢E+05 
Pe 250E+00 5.850E+05 
1.750E+00 6.171E+05 
2eeoUEFOO) 6.376E405 
Zev ouerot 6.557ER05 
Becoouetrou 6.057Et05 
eesuero0 5.715Er05 
G.250E+00 5.967E+05 
Ge, oUETOO 5.315E+05 

VARIABLE Hl 
MINVAL= 3.000E+05 
MAXVAL= 1.000E+06 
CELLAV= 5.791E+05 


+ + + + + + + + + 


K++ t+ t+etet 
t+++eete+ 


ooooo°qoeeo°oor 
NW An 
oo°o eo 
° ba os ° 
° a ° 


+ + + + ~ + + + + 

0 nal me eo 4 ne a. a 26 ny Se ae 
THE ABSCISSA IS if MIN= 2.50E-01 MAX= 4.75E+00 
HK HHH HHH HHH HHH HEHE HEHEHE HEHEHE HH HHH HHH HHH KKK HHH KH HM HK HK KKK HK KKK 


S 

oO 
4 
o++ + 


H vs Y~-Axis At LY=1 To IY=10 
For IX=1 And IZ=20 


Eeoure C22. 


1G 


PATCHCPE Oh? »PROFIL, hg Dine 15> 10> 2205 ee iy 
PLORCr lOten »W1l = , 2.000E+01, 1.750E&+02) 
TABULATION OF ser tks, AND ORDINATES... 
if v4 
2.500E-01 4.170E+01 
7.500E-01 5.000E&+01 
1.250E+00 5.790E+01 
1.750E+00 6.663£+01 
2.250E+00 7.281E+01 
2.750E+00 7.200£+01 
3.250E+00 6.108E+01 
3.750E+00 5.297E&+01 
G.250E+00 4.782E+01 
G4.750E&+00 4.620E+01 
VARIABLE bi] 
MINVAL= 2.000E+01] 
MAXVAL= 1.750&+02 
CELLAV= 5.691E+01 
i OO ee en eS ech, SLE Be aS oe aor 
0.90 + + 
0.80 + + 
0.70 + + 
0.60 + + 
0.50 + + 
0.40 + + 
0.30 + + 
0.20 + W 
0.10 W + 
0.00 Fe 0 le Fe Wow Wetts a ae oY sen a aoe gage 
0 at az aS ~4 es .6 ag, .8 9° ea 
THE ABSCiSehers MIN= 2.50E-01 MAX= 4.75£E+00 


Pagune “C7 3. 


W vs Y-Axis At IY=1 To IY=10 


For IX=1 And IZ=20 


qa 


an 
J 


1 


PATCHCPLOT8 pale te Cie los 1, 5; 5; Secu, alr, 1, 
PeOrCPLOTS ,Hl > Ss UCU E EOS mt OC CETTE) 
TABULATION OF ABSCISSA AND ORDINATES... 


xX H] 
PeoUUE=01 63762105 
PeoOCE=0) 6.165005 
Peeouety 45. 890Gr0S 
PeroUetOU 5.679E705 
@acovat00 5.5/77EP05 
VARIABLE H1 
MINVAL= 3.000E+05 
MAXVAL= 1.000E+06 
CELLAV= 5.937E+05 
LW) Sos 605 Oe een ee ere ee i 
io 0 6+ + 
u50 + cf 
tee + a 
0.60 + 1 
0.50 H—————____H__ ai 
0.40 + rs: ; eee 4 
weo0 + 1 
tec + + 
v2.10 + + 
LL 92.5 NER Se ee ere ee eS ee 
0 oat bie: ao 4 we =o eis ae ee ee 
THE ABSCISSA IS uN 2 SOE SOMA X=) 2.256700 


Figure C.24. H vs X-Axis At IX=1 TO IX=5 
For IY=5 And IZ=20 


a4 


PATCH Ur RO »PROFIL, l, 5> 5, 55 20,5.) 208 ie 
PLOT CRG? phe ec OOO TOL ley SUT Oe 
TABULATION OF ABSCISSA AND ORDINATES... 

x W1 


2 SOUESU 72s eso 
7.500E-01 6.474E+01 
lL. 2o0Ef00  S-o70ez01 
1.750E+00 4.524E+01 
2.250E+00 4.084E+01 
VARIABLE W1 
MINVAL= 2.000E+01 
MAXVAL= 1.750E+02 
CERLAV= 5. 558Er01 
L0G" 22 Ror Rea vis tig oa eae te oe ee area Se 5 ea et co 
OP as : t 
0.80 + + 
OP Uae - 
0.60 + 0 
OA Ua 2 25 
0.40 + - 
0.30 a 
ices + 
tO: : 
Cet CR en, ene rng ene nee us Cars Soe oto 0 
0 rl Te ae 4 “5 6 ne is 29 le 
THE, ABSGISSA 15 X .  MIN= 2.50E-01 MAX= 2. 25bg0 


HH KK HK HK HK HH HK HH KH HHH HEH HH KK HEH HHH KK KKK HHH HH HH HE KKK KK KKK KKKKKKK 


Figure C.25. W VS@xX-AxiceAe lon =) 
For ITY=5 And IZ=20 


1 


12 


PATCHCPLOT14 ,PROFIL, l, Soc cl, <2o, thes 
BeencPLoOnTl4 ,H1l , 5.000E+O05, 6.000E+05) 
TABULATION OF ABSCISSA AND ORDINATES... 
Z Hl 
1.Q000E+00 5.429E+05 
Cese0Eetu0 § .S35E+05 
SeUSGETU0 5. 6556705 
imo? SEU 5.895E+05 
1.845E+01 5.905E+05 
Zeeoiesie+0] 5-66 10E4+05 
VARIABLE Hl 
MINVAL= 5.000E+05 
MAXVAL= 6.Q0Q00E+05 
CELL AV=— 5.66 5E+05 
moo +. AC, Se ae SUE 2 eee ie ft. Pos aa ae + 
0.90 + + 
i260 + a 
15 Use 
Qeo0 + ne 
oo + ae 
0.490 G2 
tees 0 + + 
ieee + 7 
10 + ar 
0 ae a le ei Par Be a ee eo a 
0 Ail ee me .4 gD “6 ei 26 9s. 0 
ie ABSCISSA IS Z MIN= 1.00E+00 MAX= 2.31Et01 


HK KKK KKH HK HHH HHH HK HHH HHH HK HHH HK HK HK HHH HK HK KK KKK KK KKK KK KKKKKAKAKK 


Bigure C.26. @hevs Z-Axis At [Z=21 To I2=26 


For IX=1 And ILY=13 


see: 


1) 


PATCHCUPLONIS. 2RROFELL, l, ly 133. 13, 2iecee iy 
PLOULGR ROTTS wv! > L20D0E+ G1 34. SUCE+0 
TABULATION OF ABSCISSA AND ORDINATES... 
Z V1 
TVOUGE+GO0 1. 2e6eEs0! 
4.350E+00 2.973E+01 
9-0S0Et00 252256101 
Ies7 S5E+tO0!] 3. 062ZEr0)] 
1.8465E+01 2.919E+01 
2 OSEt0L: 2¥o2erol 
VAR EAD le Wl 
MINVAL= 1.000E+01 
MAXVAL= 4.500E+01 
CEELAV= 2.5947E+01 
1.00 +. + +....¢....+ t+. + + + + + 
0.90 + + 
[Oe 0) a + 
0.70 + + 
0.60 + V V oF 
0.50 + ere NINE IER (pre 
0.40 a ~ 
Ui sU V + 
Ors Ua: + 
0.10 V 4. 
0. a, ee eee, eee ee ne, 
0 a | a2 5 4 5 ac s7 8 9 #+1.0 
THE AbSelssclwns Z MIN= 1.00E+00 MAX= 2.31E+01 


V vs Z-Axis At IZ=21 To IZ=26 
For IX=1 And IY=13 


Figure C.27. 
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Figure C.29. V vs X-Axis At IX=1 to IX=5 
For IY=13 And IZ=23 
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